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PREFACE 

Surface  waves  have  been  found  useful  in  many  antenna  applications  and  in 
some  transmission  line  applications  for  the  past  seventy  years.  The  surface 
wave  transmission  line  is  exceptionally  efficient.  This  very  fact  makes  the 
surface  wave  a  nuisance  when  it  is  not  the  desired  mode.  Undesired  coupling 
between  components  of  a  system,  particularly  at  microwave  frequencies  and 
above,  may  occur  via  surface  waves.  Unpublished  research  by  Dr.  Gary  L. 
Johnson  into  the  attenuation  of  surface  waves  guided  by  electrically  thin 
lossy  dielectric  slabs  showed  that  the  attenuation  does  not  increase  mono- 
tonically  with  increasing  dielectric  loss.  This  thesis  extends  those  results 
to  electrically  thick  slabs.  It  is  found  that  the  same  effect  appears  for 
relatively  thick  slabs  but  not  for  extremely  thick  slabs,  an  abrupt  change 
from  decreasing  to  increasing  attenuation  occurring  at  some  thickness. 
Furthermore,  the  existence  of  an  extraordinary  mode  has  been  discovered. 

I  wish  to  express  my  gratitude  to  Dr.  Gary  L.  Johnson,  my  thesis 
advisor;  to  the  National  Aeronautics  and  Space  Administration  for  financial 
support  in  the  form  of  a  NASA  traineeship  which  made  this  research  possible; 
to  Dr.   W.  W.  Koepsel  and  Dr.  F.  W.  Harris  of  my  advisory  committee  for  their 
assistance;  to  Mr.  J.  H.  Peterson  for  his  valuable  programming  assistance; 
and  to  the  Collins  Radio  Company  and  Mr.  G.  W.  Deen  for  the  programming 
experience  which  was  invaluable  to  me  during  my  research. 
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CHAPTER  I 
INTRODUCTION 

1.1   Statement  of  the  problem.  During  the  past  seventy  years,  consid- 
erable interest  has  been  exhibited  in  surface  and  leaky  waves.  Antennas 
utilizing  surface  waves  and  leaky  waves  have  narrow  beamwidth  and  can  easily 
be  scanned  electronically;  a  large  quantity  of  literature  has  appeared  con- 
cerning the  design  and  characteristics  of  such  antennas.  A  rather  extensive 
bibliography  of  articles  published  prior  to  1961  is  given  by  Zucker  [1961]. 
There  are  cases  of  interest  in  which  surface  waves  may  be  present  but  are  not 
desirable.  For  example,  radio  frequency  interference  (RFI)  is  propagated 
along  electric  power  lines  by  a  surface  wave;  at  microwave  frequencies  and 
above,  the  use  of  dielectric  substrates  with  dimensions  comparable  to  a  wave- 
length makes  the  excitation  of  surface  waves  \/ery   likely.  It  appears  that 
the  introduction  of  losses  should  produce  attenuation  and  thus  help  to 
eliminate  such  waves.  Little  information  appears,  however,  concerning  the 
lossy  case  for  a  thick  dielectric  coating.  Goubau  [1950]  and  Collin  [1960] 
treat  low  loss  cylindrical  conductor  cases;  Barlow  and  Cullen  [1953]  and 
Johnson  [1966]  treat  the  thin  planar  slab  lossy  case. 

It  is  the  purpose  of  this  research  to  present  data  from  a  computer  aided 
numerical  analysis  to  demonstrate  the  effect  of  lossy  materials  upon  the 
propagation  of  surface  waves  over  thick  dielectric  slabs.  An  ideal  planar 
case  is  treated  in  which  it  is  assumed  that  the  structure  considered  is 
infinitely  wide,  that  the  source  and  other  discontinuities  are  sufficiently 
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removed  from  the  region  under  consideration  that  no  reflected  waves  are 
present,  and  that  the  system  is  in  the  steady  state,  the  source  having  an 
exp(ju)t)  time  dependence.  Only  TM  modes  having  no  field  variation  along 
the  width  of  the  structure  are  considered.  The  materials  are  allowed  to  be 
lossy,  but  are  assumed  to  be  linear,  homogeneous,  isotropic,  and  time- 
invariant. 

The  dielectric  slab  is  coated  onto  the  surface  of  a  perfect  conductor. 
It  can  be  shown  that  the  use  of  a  good  but  not  perfect  conductor  has  a  neg- 
ligible effect.  If  it  is  desired  to  consider  the  conductor  losses  a 
perturbational  analysis  should  be  entirely  satisfactory.  The  assumption 
that  there  is  no  reflection  in  the  lengthwise  or  axial  direction  is  taken 
merely  for  convenience  since  the  case  having  reflection  can  be  handled  with 
the  reflectionless  fields  and  the  usual  transmission-reflection  techniques. 
The  least  realistic  assumption  is  likely  to  be  the  assumption  that  the  slab 
is  of  infinite  width,  which  is  equivalent  to  the  assumption  of  no  reflection 
from  discontinuities  at  the  edges  of  a  finite  slab  or  that  the  region  under 
consideration  is  sufficiently  removed  from  the  discontinuities  that  losses 
eliminate  the  reflected  wave.  Techniques  for  reducing  the  reflection  at  the 
edges  of  a  finite  slab  are  described  by  Zucker  [1961].  The  consideration  of 
a  finite  slab  width  would  make  the  problem  a  rather  formidable  boundary  value 
problem. 

1.2   Discussion  of  the  terms  "Surface  Wave"  and  "Leaky  Wave."  It  is 
convenient  to  first  describe  surface  and  leaky  waves  for  the  lossless  case. 
At  a  boundary  between  two  different  media  there  must  result  reflection  and 
transmission  of  any  electromagnetic  wave  which  is  present.  The  surface  wave 
in  this  lossless  case  is  totally  reflected  at  the  boundary  so  that  the  field 
in  the  less  dense  medium  is  an  evanescent  field  and  the  field  in  the  more 


dense  medium  exhibits  a  standing  wave  behavior.  The  direction  of  propaga- 
tion is  parallel  to  the  interface.  There  is  no  energy  propagating  into  the 
less  dense  medium  from  the  more  dense  medium  or  vice-versa.  The  wave  is 
unattenuated  in  the  direction  of  propagation  but  decays  exponentially  with 
increasing  distance  from  the  interface  in  the  less  dense  medium. 

The  leaky  wave  is  one  having  no  real  critical  angle  so  that  the  wave 
cannot  be  totally  reflected  at  the  interface.  Instead  energy  now  passes 
from  the  more  dense  medium  to  the  less  dense  medium,  and  as  a  result  the 
amplitude  of  the  wave  decreases  along  the  interface.  The  wave  in  the  less 
dense  medium  is  still  an  evanescent  wave  although  the  direction  of  propaga- 
tion is  tilted  into  the  less  dense  medium.  The  amplitude  of  the  wave 
increases  without  bound  as  the  distance  into  the  less  dense  medium  increases; 
though  this  is  not  a  physically  meaningful  wave  in  this  ideal  case,  the 
increasing  amplitude  is  caused  by  an  energy  source  infinitely  removed  from 
the  region  under  consideration.  Such  a  source  is  not  present  in  physically 
real  situations;  the  resulting  fields  from  the  physically  real  source  are 
known  to  be  described  reasonably  well  by  the  ideal  model  at  reasonable 
distances  from  the  interface. 

Qualitatively  speaking,  the  presence  of  losses  modifies  the  direction  of 
propagation  in  such  a  direction  as  to  support  the  losses.  If  the  less  dense 
region  is  lossy  the  wave  in  that  region  will  not  be  evanescent.  The  wave 
will  not  be  totally  reflected  at  the  interface  except  possibly  for  certain 
relations  between  the  complex  permittivities  and  permeabilities  of  the  two 
media.  In  any  lossy  case,  both  the  surface  wave  and  the  leaky  wave  must 
diminish  in  amplitude  as  the  wave  propagates  along  the  interface.  The 
distinction  between  a  surface  wave  and  a  leaky  wave  in  a  lossy  case  is  in 
the  direction  of  attenuation. 


CHAPTER  II 

THE  CHARACTERISTIC  EQUATION 

2.1   Derivation  of  Equations.  Consider  the  geometry  of  Figure  1, 
showing  the  dielectric  slab  and  the  perfect  conductor  to  which  the  slab  is 
attached,  and  the  orientation  of  the  cartesian  coordinate  axes.  The  y-axis 
is  directed  out  of  the  page;  the  slab  and  conductor  are  infinite  in  extent 
in  the  y-direction  and  in  the  z-direction.  The  slab  thickness  is  t,  with 
the  conductor-slab  interface  in  the  y-z  plane.  Both  the  dielectric  slab 
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Figure  1.  Geometry  of  the  Problem. 


and  the  surrounding  medium  are  linear,  homogeneous,  isotropic  and  time- 
invariant;  thus  y, ,  y2,  e-i ,  and  e2  are  scalar  complex  constants.  It  is 
assumed  that  y^eo  *  ^\z\   anc'  e?  >  eT 


The  losses  appear  in  the  permittivity  and  permeability  with  the 
following  notation: 


£  =  ere0  =  e0(er'  "  jerM) 


vniv'   -  jy  ") 


rM0   M0VMr 

In  these  expressions  e  and  y  are  the  permittivity  and  permeability  of  the 
material  while  the  zero  subscripts  denote  the  free  space  values,  e  '  is 
the  relative  AC  capacitivity  or  dielectric  constant  while  e  "  is  the  rela- 
tive dielectric  loss  factor.  Similarly  y  '  is  the  relative  AC  inductivity 
and  y  "  is  the  relative  magnetic  loss  factor. 

In  writing  the  field  quantities  the  conventions  employed  by  Harrington 
[1961]  will  be  used.  The  time  variation  e^w  is  suppressed.  Since  only 
the  sinusoidal  steady  state  is  considered,  there  is  no  requirement  to  con- 
sider the  actual  instantaneous  values.  As  used  in  this  thesis  the  vector 
magnetic  potential  A  will  be  defined  by  the  relations 

ft  =  v  x  A  [2.1.1] 

-jcoeO  =  v  •  t  [2.1  .2] 

where  $  is  a  scalar  electrical  potential.  This  choice  causes  A"  to  satisfy 
the  well-known  Helmholtz  equation.  Subscripts  1  and  2  will  be  used  to  iden- 
tify the  field  quantities  in  the  surrounding  medium  and  in  the  dielectric 
slab  respectively. 

Before  beginning  the  derivation  of  the  characteristic  equation,  an  out- 
line of  the  procedure  used  will  be  presented.  An  elementary  plane  wave 
function  for  each  of  the  two  homogeneous  regions  0  <_  x  <_  t,  t  <_  x  will  be 
assumed.  These  functions  will  be  used  for  the  vector  magnetic  potential. 


The  resulting  field  components  will  be  calculated  from  Maxwell's  equations, 
and  boundary  conditions  will  be  matched  at  x  =  t.  Finally,  two  additional 
necessary  relations  will  be  obtained  from  the  separation  equation. 

It  is  assumed  that  the  vector  magnetic  potential  is  given  by  the 
expressions 

A2  =  C2  sin(ux)  e    az    O^x^t  [2.1 .3] 

-vx  -jk  z 
A1  =  C]  e   e  L   az        t  <  x  [2.1.4] 

t     =  0  x  <_  0  [2.1.5] 

where  C,  and  C2  are  undetermined  coefficients.  Thus  in  the  dielectric  the 
x-component  of  the  propagation  constant  is  u,  while  in  the  surrounding 
medium  it  is  jv.  The  z-component  of  the  propagation  constant  is  k  in  both 
regions  as  is  required  by  Snell's  Law.  The  fields  are  zero  beneath  the 
surface  of  the  perfect  conductor;  it  is  immaterial  how  thick  this  conductor 
is  or  what  lies  below  it,  as  any  -x  traveling  wave  is  totally  reflected  at 
the  surface  in  this  ideal  case.  H  is  calculated  from  equation  2.1.1: 

-jk  z^ 
ft2  =  vx^2  =  -uC2  cos(ux)  e    a  [2.1.6] 

-vx  -jk  z 

^  =  vx^  =  v^e   e   z  a  [2.1.7] 

E  is  calculated  from  Maxwell's  equations  for  this  source-free  case: 

l  uk7  "Jk7Z->  n2  -J'k7z+ 

E9  =  —  vxH9  = ^C0cos(ux)e  a     +  -^ — C0sin(ux)e  a,  [2.1.8] 

2       y9         2  we0     2        K      '  x       jioe9     2  z 


,  vk  -vx  -jk  z  w2  -vx  -jk  z 

f-l-vxH._i.re      e      zl    -    v      c  e      zt  [2J>g] 

I    y-,  I    we,   I  X    Jooe-i   I  Z. 

At  x  =  t,  the  tangential  components  of  both  £  and  ft  must  be  continuous;  con- 
tinuity of  E  requires  that 

C  C 

—  u2  sin(ut)  =  -  —  v2e~vt  [2.1.10] 

£2  el 


while  continuity  of  H  requires  that 

-C2u  cos(ut)  =  C,ve"vt.  [2.1.11] 


The  continuity  of  normal  electric  flux  leads  to  2.1.11  also,  so  it  imposes  no 
additional  constraints.  The  necessary  boundary  condition  of  zero  tangential 
component  of  E  at  the  conductor  surface  is  satisfied  by  E  ?.  Taking  the 
ratio  of  2.1.10  to  2.1.11  to  eliminate  C,  and  C2,  we  have 


^-tan(ut)  =  — ^  ; 
e2         el 


it  is  convenient  to  rearrange  this  to  yield 


vt  =  —  ut  tan(ut).  [2.1.12] 

er 


In  equation  2.1.12  the  relative  permittivity  between  the  two  media  has  been 
defined  as 

£r  =  e2/eT 
Equation  2.1.12  is  the  characteristic  equation  relating  u  and  v. 


Another  equation  relating  u  and  v  is  necessary  in  order  to  obtain  a 

solution.  The  Helmholtz  equation  is  a  direct  consequence  of  Maxwell's 

equations,  and  it  has  been  separated  in  the  cartesian  coordinate  system, 
yielding  the  separation  equation 

kx  +  ky  +  kz  =  k2  [2.1.13] 

where  k  is  the  propagation  constant  or  wave  number  of  the  medium  and  k  ,  k  , 
and  k  are  the  components  of  the  vector  propagation  constant  of  the  wave.  In 
the  case  considered  here,  k  is  zero  in  both  regions,  k  is  the  z-component 
in  both  regions,  jv  is  the  x-component  in  the  surrounding  medium,  and  u  is 
the  x-component  in  the  dielectric.  Therefore  the  application  of  2.1.13  to 
each  medium  yields  the  two  equations 

u2  +  k2  =  k2  =  u2v2e2  [2.1.14] 

(jv)2  +  k2  =  k2  =  u^e-,  [2.1.15] 

containing  the  three  unknowns  u,  v,  and  k  .  k  can  be  eliminated  by  subtrac- 
tion of  2.1.15  from  2.1.14  to  yield 

U2  +  V2  =  u)2(y2£2  -  y^-,).  [2.1.16] 

It  is  convenient  to  re-write  this  as 

(ut)2  +  (vt)2  =  (2tt[^]  SVr2er2   -  Mrlerl)2  [2.1.17] 

where  Ao  is  the  free-space  wavelength  of  the  wave,  and  u  , ,  u  „,  e  , ,  z  2   are 
the  various  permeabilities  and  permittivities  normalized  to  the  free  space 


values.  The  unknown  v  may  be  eliminated  from  equations  2.1.12  and  2.1.17  to 
yield 

(ut)2  +  1-2-  (ut)2  tan^ut)  =   (2Ti[^]/Mr2cr2  -  yrl^rl)2.  [2.1.18] 

r 

This  is  a  non-linear,  transcendental  equation  with  coefficients  which  are  in 
general  complex.  This  is  the  equation  which  will  be  solved  to  obtain  u. 

2.2   Graphical  Analysis  of  Real  Roots  in  the  Lossless  Case.  For  the 
lossless  case,  the  parameters  of  the  surface  wave  modes  can  be  obtained  from 
a  graphical  solution.  Since  e  ,,  e  «i  y  i,  and  y  ?  are  real  for  the  lossless 
case,  and  k  and  v  must  be  real  for  a  surface  wave  mode  in  the  lossless  case, 
u  is  also  real.  This  is  shown  by  considering  Equations  2.1.14  and  2.1.15. 
When  u  and  v  are  real  they  represent  only  two  unknowns  whereas  they  would 
represent  four  unknowns  in  the  complex  case.  Two  simultaneous  equations  in 
two  unknowns  can  be  solved  by  graphing  the  two  relationships  and  finding  the 
points  of  intersection  as  is  commonly  known.  Equations  2.1.12  and  2.1.17 
are  well -suited  for  this  purpose;  Equation  2.1.17  is  recognized  as  a  circle 
in  the  ut-vt  plane,  while  Equation  2.1.12  has  the  general  features  of  a 
tangent  function. 

Figure  2  demonstrates  the  graphical  solution  for  e  «  =  2.0,  e  ■,   =1.0, 
y  2  =  v  -i   -   1.0,  and  several  values  of  t/Ao.  The  radius  of  the  circle  is 
2(t/Xo)TT  for  this  set  of  parameters.  Modes  (solutions  to  2.1.18)  are 
numbered  in  order  of  increasing  ut  for  the  lossless  case.  The  same  mode 
number  is  retained  as  losses  are  added.  The  TMo  mode  refers  to  the  real 
intersection  of  the  circle  and  the  quasi-tangent  function  in  the  interval 
0  <  ut  <  tt/2.  This  mode  always  propagates  as  a  surface  wave  regardless  of 
how  thin  the  slab  is.  Higher  order  modes  may  represent  surface  waves,  leaky 
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Figure  2.  Graphical  Solution  for  the  Lossless  Case, 
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waves,  or  non  physical  waves  depending  upon  the  radius  of  the  circle  and 
losses. 

It  is  clear  that  there  is  either  one  or  no  real  intersection  of  the  two 
curves  between  any  two  successive  integral  multiples  of  tt/2  on  the  ut  axis. 
In  general,  there  will  be  exactly  one  real  intersection,  and  hence  one  sur- 
face wave  mode,  having  a  value  of  ut  between  I<tt/2  and  ( k+1  )tt/2 ,  where  k,  the 
mode  number,  is  an  even  integer,  provided  that  the  circle  radius  is  greater 
than  or  equal  to  k-n/2.     All  higher  ordered  surface  wave  modes  have  a  defi- 
nite cutoff  frequency  below  which  they  become  leaky  wave  modes.  Cutoff  for 
the  k   mode  occurs  when  the  circle  radius  is  equal  to  k-rr/2. 

Obviously  there  may  also  be  intersections  in  the  intervals 
m-n/2  <_  ut  <_  (m+1  )tt/2  where  m  is  an  odd  integer  if  the  radius  of  the  circle 
is  at  least  miT/2.  There  may  be  either  two  distinct  intersections  or  a  point 
of  tangency  in  an  interval  in  which  the  radius  of  the  circle  lies;  for  the 
intervals  such  that  rmr/2  <_  ut  <  (m+l)iT/2  and  the  circle  radius  is  at  least 
(m+1 )%/2   there  will  be  exactly  one  distinct  intersection.  The  value  of  v 
which  results  in  these  intervals  is  negative,  so  that  the  fields  would 
increase  in  magnitude  with  increasing  distance  from  the  slab.  There  would 
not  be  any  radiation  from  the  slab.  This  describes  neither  a  surface  wave 
nor  a  leaky  wave,  so  such  modes  will  be  regarded  as  non-physical. 

Negative  values  of  ut  are  not  shown  in  Figure  2.   It  is  apparent  from 
Equations  2.1.12  through  2.1.18  that  if  u-,  is  a  solution,  -u-,  is  also  a  solu- 
tion resulting  in  the  exact  same  values  of  v  and  k  .   Inspection  of  Equations 
2.1.6  through  2.1.11  shows  that  if  u  changes  sign  but  v  remains  unchanged 
then  either  C-,  or  C~  must  change  sign.   If  C~  changes  sign  there  is  no  dif- 
ference between  the  fields  resulting  from  u,  and  -u-, .   If  C-,  changes  sign, 
all  field  components  are  reversed  in  direction;  this  corresponds  to  the 
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original  field  pattern  viewed  half  a  guide  wavelength  in  either  direction 
from  the  original  point.  Therefore,  the  consideration  of  negative  values  of 
u  adds  no  new  information  and  is  unnecessary. 

Considerable  insight  into  the  behavior  of  the  propagation  constants  can 
be  obtained  with  this  graphical  technique.  In  order  to  determine  the  behav- 


ior of  k  ,  Equations  2.1.14  and  2.1.15  are  also  necessary.  These  two  equa- 
tions are  easily  rearranged  to  show  that  k-.  <_  k  <  k?  for  the  cases  discuss 
in  this  section,  that  is,  cases  in  which  u,  v,  and  k  are  real. 


ed 


Consider  cases  in  which  Ao,  the  free  space  wavelength,  is  the  only 
parameter  allowed  to  vary  in  Figure  2.  As  Ao  is  changed,  only  the  radius  of 
the  circle  changes.  As  Ao  decreases,  t/Ao  increases  and  ut  and  vt  increase. 
Therefore,  u,  v,  and  k  all  increase  since  t  is  fixed.  For  ^ery   large  Ao, 
so  that  the  radius  is  small  compared  to  tt/2,  ut  is  approximately  equal  to 
the  radius  and  vt  is  small,  k  tends  toward  its  lower  bound  k, .  Relatively 
speaking,  the  wave  will  be  loosely  bound,  and  the  wavelength  close  to  A-,. 
For  yery   small  Ao,  so  that  the  circle  radius  is  much  greater  than  tt/2,  ut  is 
approximately  tt/2  for  the  first  mode,  and  the  corresponding  vt  is  approxi- 
mately the  radius  of  the  circle,  k  tends  toward  its  upper  bound  k?.  Again, 
relatively  speaking,  the  wave  will  be  tightly  bound,  and  the  wavelength  close 
to  Ao- 

The  behavior  of  any  higher  propagating  modes  is  obviously  similar,  ut 
being  approximately  the  radius  and  vt  being  small  for  a  case  slightly  away 
from  the  cutoff  region,  ut  being  near  (slightly  less  than)  an  odd  multiple 
of  tt/2  and  vt  being  somewhat  less  than  the  radius  for  a  very  large  value  of 
radius. 

2.3   Some  Notes  Concerning  Solutions  in  More  General  Cases.  No  closed 
form  solution  of  Equation  2.1.18  is  known.  The  graphical  technique  of 
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section  2.2  is  not  useful  when  the  solutions  are  complex  since  there  are 
four  unknowns  rather  than  two.  An  iterative  numerical  procedure  appears  to 
be  the  most  straightforward  method  of  solution.  In  order  to  use  this  tech- 
nique it  is  necessary  to  have  a  reasonable  estimate  of  the  solution  to  use 
as  a  starting  value.  In  particular,  it  would  be  appropriate  to  be  able  to 
place  bounds  upon  the  value  of  u  for  a  given  mode. 

As  a  first  step  in  examining  the  complex  case  it  is  convenient  to 
introduce  the  substitution  z  =  ut;  note  that  this  variable  z  is  not  the 
coordinate  value  z,  but  another  distinct  variable.  Equation  2.1.18  then 
becomes 

Z2(1    +  I_  tan*(z))   =    (27r[^]    /yr2cr2   -   Mrlsrl)2  [2.3.1] 

or 


yr2£r2  "  yrl£rl       = 


z  -  27Ter[y   /    e/  +  tan^(z)  "  u'  [2.3.2] 


If  f, (z)   is  defined  as 


t         /yr2£r2  "  yr1£rl 
f^z)  =  z  -  2TT£r[y   /    Ep2  +  tan^(z)  [2.3.3] 

the  solutions  desired  are  the  zeros  of  f-.(z).     For  convenient  reference, 
equations  2.1.12,  2.1.14,  and  2.1.15  are  also  repeated  here: 

vt  =  1_  z  tan(z)  [2.3.4] 

er 

u2  +  k2z  =  k|  =  u2v2z2  [2.3.5] 

-v2  +  k2  =  k2  =  o)2yl£l.  [2.3.6] 
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In  obtaining  equation  2.3.2  from  equation  2.3.1  at  least  one  root  has 
been  ignored  as  a  consequence  of  the  convention  that  the  radical  sign  indi- 
cates the  root  with  the  non-negative  real  part,  or  the  root  with  the  posi- 
tive imaginary  part  if  both  roots  are  purely  imaginary.  This  does  not 
detract  from  the  generality  of  the  resulting  solution  as  has  been  discussed 
in  section  2.2;  i.e.,  while  solutions  with  negative  real  parts  are  distinct 
mathematical  solutions,  they  are  not  distinct  physical  solutions.  Having 
restricted  the  choice  of  z  (and  hence  u)  a  single  unique  value  of  v  is  thus 
determined  by  equation  2.3.4;  therefore  one  of  the  values  which  results  from 
equation  2.3.6  is  extraneous,  k  ,  however,  is  not  so  restricted  mathemati- 
cally; the  two  values  k  =  +/k£  -  uz  are  both  valid  physically.  Clearly, 
the  significance  of  the  existence  of  two  values  is  that  propagation  may 
occur  in  either  the  +z  or  the  -z  direction,  or  both  (when  there  is  reflec- 
tion, say).  Since  the  features  of  the  two  waves  are  identical  except  for 
the  direction  of  propagation  it  will  be  sufficient  to  examine  only  the  +z 
propagating  solution  in  this  research;  if  the  more  general  case  must  be 
considered  the  appropriate  change  of  signs  may  be  used  to  describe  the 
oppositely  traveling  wave. 

Consider  now  the  meaning  of  the  signs  of  the  three  propagation  constants 
u,  v,  and  k  .  As  the  wave  inside  the  dielectric  slab  is  a  standing  wave,  the 
signs  of  the  real  and  imaginary  parts  of  u  (and  z  =  ut)  do  not  have  any 
physical  significance.  The  signs  of  the  real  and  imaginary  parts  of  v  and  k 
determine  the  direction  of  propagation  and  the  direction  of  attenuation  as 
indicated  in  table  1.   For  a  surface  wave  traveling  in  the  +z  direction,  it 
is  required  that  k  lie  in  the  fourth  quadrant  of  the  complex  k  plane,  and 
that  v  lie  in  the  right  half  of  the  complex  v  plane.  For  a  leaky  wave  mode 
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k  must  still  lie  in  the  fourth  quadrant,  but  v  must  lie  in  the  left  half  of 


the  complex  v  plane. 


phase  travel 

Re[kz] 

in 

+z 

direction 

attenuation 

Im[kz] 

in 

-z 

direction 

attenuation 

Re[v] 

in 

+x 

direction 

phase  travel 

Im[v] 

in 

+x 

direction 

TABLE  1 
SIGNIFICANCE  OF  ALGEBRAIC  SIGNS  OF  v  AND  k  . 

positive  negative 

in  -z  direction 
in  +z  di  recti  on 

in  -x  direction 
in  -x  direction 

It  is  now  informative  to  consider  the  vector  propagation  constant.  As 
defined  by  Harrington  [1961]  the  vector  propagation  constant  y  =  a  +  jB 
represents  a  wave  having  equi phase  surfaces  normal  to  3  and  equi amplitude 
surfaces  normal  to  a.  In  this  research  a   and  6  in  region  1  are  given  by  the 
expressions 


a  =  v'ax-kz"az  [2.3.7] 


=  v"ax+kz'az.  [2.3.8] 


Since  these  vectors  are  independent  of  the  coordinates  (x,y,z)  the  wave  in 
region  1  is  a  plane  wave  traveling  in  the  3  direction  and  attenuating  in  the 
a  direction.  The  wave  in  region  2  is  a  standing  wave  composed  of  two  plane 
waves  which  are  totally  reflected  at  the  perfect  conductor  and  are  reflected 
also  at  the  interface  between  the  two  media.  For  a  lossless  surface  wave 
total  reflection  occurs  at  the  dielectric  interface  so  that  the  wave  in 
region  1  is  an  evanescent  wave.  Examination  of  the  behavior  of  a  and  6  is 
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of  some  interest,  particularly  their  directions.  The  angle  eo  between  the 

B 

B  and  a„  directions  and  the  angle  e  between  the  J  and  a  directions  are 

z  3    a  z 

given  by  the  expressions 


e6  =  tan_1(v7kz')  [2.3.9] 

e  =  tan"1(v7(-k'1)),  [2.3.10] 

a  z 

where  the  notation  tan-  denotes  the  arctangent  function.  Figure  3  shows  a, 
B»  e  ,  and  e  for  surface  waves  over  lossy  slabs  and  leaky  waves.  As 
mentioned  earlier,  the  surface  wave  propagates  into  the  lossy  dielectric  in 
the  ?  direction  to  account  for  the  losses.  Since  a  and  B  are  orthogonal  for 
region  1  lossless,  and  a  must  be  in  or  on  the  boundary  of  quadrant  1,  then 
B  will  be  in  quadrant  4.  A  similar  analysis  for  leaky  waves  shows  that  a 
must  lie  in  the  fourth  quadrant  and  B  in  the  first  quadrant. 

The  surface  wave  modes  of  the  lossless  case,  i.e.,  the  real  solutions, 
have  already  been  considered.  The  leaky  wave  modes  of  the  lossless  case 
will  now  be  considered.   If,  for  a  given  even  mode  number,  there  is  no  real 
intersection  of  the  circle  and  the  quasi-tangent  function,  the  corresponding 
root  of  2.1.18  will  be  complex  and  describe  a  leaky  wave  mode.  A  straight- 
forward examination  of  equations  2.3.5  and  2.3.7-10  shows  for  this  case, 
that  u  and  z  lie  in  the  first  quadrant,  k  lies  in  the  fourth  quadrant,  and 
v  lies  in  the  second  quadrant  (of  their  respective  complex  planes).  Although 
the  restriction  of  u  to  the  first  quadrant  is  informative,  an  entire  quadrant 
is  still  a  rather  large  region  to  scan  for  a  particular  root.  It  seems 
intuitively  obvious  that  a  given  leaky  wave  mode  solution  must  lie  in  a  rel- 
atively small  region  as  in  the  case  of  a  given  surface  wave  mode;  however, 
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Surface  Wave  Over  Lossy  Dielectric 


Leaky  Wave 


Figure  3.  Typical  Directions  of  Components  of  the  Vector 
Propagation  Constant. 
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no  mathematical  description  of  such  a  region  has  been  determined  by  the 
author.  This  region  is  obviously  a  subset  of  the  set  [z |Re(ztan(z))<0]. 

Similarly,  it  is  extremely  difficult  to  define  the  permissible  region 
for  the  lossy  surface  wave  case  and  the  lossy  leaky  wave  case,  as  the  com- 
plex relative  permeability  between  the  two  materials  also  enters  into  the 
determination  of  the  sign  of  Re[v]. 

It  is  interesting  to  note  that,  in  the  lossless  leaky  wave  case,  if  z, 
v,  and  k  satisfy  the  equations  2.3.3  to  2.3.6,  then  their  conjugate  values 
z*,  v*,  and  k*  also  satisfy  those  equations,  though  the  conjugate  values 
describe  a  wave  which  is  not  physically  meaningful.  Using  the  first  term  of 
the  Taylor  series  expansion  of  tan(z)  in  Equation  2.3.1,  it  is  apparent  that 
the  lossless  surface  wave  mode  solutions  are  at  least  double  roots,  and  one 
might  expect  the  two  roots  to  diverge  and  proceed  along  approximately  conju- 
gate paths  as  losses  are  added. 

Another  interesting  and  relevant  matter  is  the  behavior  of  tan(z)  for 
arguments  even  moderately  removed  from  the  real  axis.  For  b=TT,  tanh(b)  - 
0.99595;  for  b=2-rr ,  tanh(b)  =  0.99999;  tanh(b)  ■+  1  monotonically  as  b 
increases  from  0.  Also,  tanh(b)  is  an  odd  function.  For  z=a+jb, 


.   /  x    tan(a)  +  jtanh(b) 
tarUzj  "  1  -  jtan(a)tanh(b) 


(a)tanh(b) 

so  that  for  tanh(b)  -   +1,  tan(z)  -   +jl.  Under  these  conditions,  f-.(z)  takes 
on  the  value 


yr2er2  "  yrlerl 


f^z)  -  Z  -  2Trer[t/Xo]  /     e/  -  1 

For  the  lossless  case,  then,  the  roots  must  clearly  be  "close"  to  the  real 
axis,  as  |f,(z)|>c  -  |Im[z]|  for  z  sufficiently  removed  from  the  real  axis 
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This  is  significant  information  for  use  with  a  numerical  technique,  as  the 
region  of  search  can  be  restricted  to  the  vicinity  of  the  real  axis. 


CHAPTER  III 

THE  COMPUTER  PROGRAM  AND  NUMERICAL  TECHNIQUES 

3.1       Development  of  the  Program.     The  early  efforts  employed  in  this 
research  utilized  various  rearrangements  of  Equation  2.1.18.     Some  of  these 
mathematically  equivalent  rearrangements  are: 


z2  +  -ij-  Z2  tan*   (z)  -   (2ir[t/Ao]  /ur2er2  -  yrleH)2     =0;  [3.1.1] 


-   /4TT2[t/Ao]2   (e-y.  -  e.y,)   -  fV  tan2   (z)         =  0; 


r2Mr2         rlMrr       JJ 


[3.1.2] 


z  tan  (z)  -  er   /^[t/Ao]2  {v^r2   -  yrlerl)  -  z2    =0.       [3.1.3] 

Unfortunately,  such  rearrangements  are  not  equally  useful  with  numerical 
techniques;  whether  or  not  a  given  iterative  procedure  converges  to  a  valid 
result  depends  upon  the  expression  iterated  upon  as  well  as  the  starting 
guess  provided.  The  mathematical  determination  of  whether  or  not  a  given 
expression  is  usuable  in,  say,  linear  iteration,  is  a  problem  of  the  same 
magnitude  as  the  analytical  solution  of  the  original  problem.   It  has  been 
found  to  be  more  practical  to  experimentally  rather  than  analytically  test 
whether  or  not  an  expression  is  satisfactory. 

Newton's  method  was  the  first  technique  attempted.   In  the  earliest 
version  the  complex  quantity  was  considered  as  a  single  variable,  even 
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though  it  was  realized  that  Newton's  method  is  not  generally  valid  for  com- 
plex variables.  Later  efforts  considered  the  problem  in  terms  of  two  real 
variables.  The  algorithms  used  did  not  yield  convergent  solutions,  so 
Newton's  method  was  discarded. 

The  second  technique  attempted  was  linear  iteration;  the  most  satis- 
factory expression  found  for  the  iteration  is 


-,   -  7       rt/i  1  /yr2£r2  "  M£rl  .  n  ,  d1 

z  -  2TT£r[t/Ao]  /  e^2  +  tan*  (z)  [3J-4] 


The  basic  problem  with  this  technique  is  that  it  is  not  self-starting,  i.e., 
the  starting  value  must  be  quite  close  to  the  root  or  the  procedure  diverges. 
Also,  it  appears  that  the  starting  value  must  be  in  a  certain  direction  in 
the  complex  z  plane  from  the  root,  and  this  direction  depends  on  all  of  the 
parameters  of  the  case  under  consideration.  Linear  iteration  was  not  com- 
pletely abandoned,  but  temporarily  put  aside  to  allow  investigation  of  other 
techniques. 

A  simple-minded  form  of  steepest  descent  was  next  tried.  In  this 
approach,  termed  the  "garbage  can  lid"  approach,  the  error,  or  magnitude  of 
the  left  hand  side  of  an  expression  f(z)  =  0,  is  calculated  at  24  points  on 
a  circle  about  some  starting  point  z.  The  new  point  is  the  point  on  this 
circle  having  the  smallest  error  associated  with  it.  Once  a  minimum  has 
been  isolated,  the  circle  radius  is  shortened  and  the  procedure  repeated 
until  the  minimum  is  located  accurately.  Hopefully,  the  minimum  will  be 
zero.  This  technique  is  yery   simple-minded,  yery   reliable,  and  requires  a 
relatively  large  number  of  iterations.  It  is  not  true  steepest  descent,  and 
the  ultimate  accuracy  is  more  limited  by  the  number  of  points  taken  on  the 
circle  than  by  machine  precision. 
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Steepest  descent  is,  in  general,  a  very   attractive  procedure,  at  least 
in  theory.  The  method  is  widely  known  and  documented  in  many  texts,  such  as 
Kunz  [1957].  A  factor  which  is  strongly  in  its  favor  is  that  it  is  simple 
to  visualize.  The  unknowns  must  be  real,  as  must  the  expression  to  be  mini- 
mized.  If  the  objective  is  the  location  of  zeros,  the  error  expression  must 
be  positive  semidefinite;  this  requirement  is  easily  met  by  using  the  magni- 
tude or  the  square  of  the  expression.  As  steepest  descent  is  a  minimizing 
procedure,  extraneous  results  can  be  obtained  if  non-zero  minima  exist,  and, 
possibly,  if  saddle  points  exist. 

For  true  steepest  descent  it  is  necessary  to  proceed  in  the  direction 
opposite  the  gradient;  it  is  necessary,  then,  to  calculate  the  partial 
derivatives  of  the  error  with  respect  to  the  unknowns,  which  unknowns  are  in 
this  case  the  real  and  imaginary  parts  of  z.   In  this  particular  case,  the 
generation  of  analytic  expressions  for  the  derivatives  is  almost  out  of  the 
question.  The  hand  calculation  of  them  is  so  tedious  and  lengthy  as  to  be 
virtually  impossible.  An  attempt  to  generate  such  expressions  using  PL/1 
FORMAC  resulted  in  expending  two  hours  of  System  360  computer  time  and 
obtaining  expressions  for  both  derivatives,  each  of  which  required  approxi- 
mately 350  cards.  Therefore  it  was  necessary  to  use  central  difference 
approximations  for  the  derivatives.  It  is  intuitively  obvious  that,  provided 
the  signs  of  the  derivatives  are  calculated  correctly,  the  procedure  will 
converge,  although  the  convergence  rate  may  be  \/ery   slow  if  there  is  a  large 
error  in  the  magnitude  of  the  derivative  approximations.  Far  from  the  root, 
or  the  minimum,  the  derivatives  may  be  determined  reasonably  accurately. 
However,  if  the  error  is  not  smooth  and  well  behaved  in  the  vicinity  of  the 
minimum,  the  numerical  calculation  may  result  in  an  incorrect  sign  as 


23 


indicated  in  Figure  4.  When  finite  difference  derivative  approximations  are 
used,  steps  must  be  taken  to  deal  with  such  cases  as  pictured  in  Figure  4, 
and  with  the  possibility  that  the  approximations  may  result  in  zero  values 
when  the  gradient  is  not  truly  zero. 


f(x) 


Central  Difference 
Approximation 


Xo-AX 


Figure  4, 


Illustration  of  the  Possibility  of  an  Erroneous 
Derivative  Approximation  Sign. 


3.2   Final  Version  of  the  Basic  Program.  Linear  iteration,  steepest 
descent,  and  the  garbage  can  lid  procedure  are  all  used  in  the  final  version 
of  the  program.  The  control  section  of  the  program,  which  performs  such 
tasks  as  reading  the  input,  choosing  the  starting  guess  for  the  unknown,  and 
writing  the  output,  is  described  in  Section  3.3. 

The  basic  program,  an  iteration  loop,  determines  the  value  of  the 
unknown  Z,  which  is  U*T  as  defined  in  Chapter  2.  A  flowchart  of  this  iter- 
ation loop  is  included  in  Appendix  A.  The  DO  statement  is  indexed  to  allow 
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a  maximum  of  100  iterations  through  this  main  loop.  Using  linear  iteration, 
a  value  ZCK  (the  left  hand  side  of  3.1.4)  is  obtained  from  the  present  value 
of  Z  (the  starting  value  on  the  first  pass),  used  in  the  right  hand  side  of 
3.1.4.  The  error  expression  is  evaluated  from  Z  and  ZCK  to  obtain  FERR  and 
FCKERR  respectively.   If  ZCK  results  in  a  smaller  error  than  Z  and  ZCK  is 
suitably  close  to  Z  so  that  the  same  root  is  being  approached,  Z  and  ZCK  are 
interchanged  so  that  the  better  approximation  to  the  root  is  used  in  the 
remainder  of  the  pass. 

The  error  is  tested  next;  if  it  is  small  enough  that  the  true  root  has 
been  satisfactorily  approximated,  the  iteration  loop  is  exited.  Otherwise, 
the  central  difference  approximations  to  the  partial  derivatives,  FPX  and 
FPY,  are  calculated,  the  gradient  DELF  is  constructed,  and  the  magnitude  of 
the  gradient,  DELMAG,  is  tested  to  be  certain  it  is  non-zero.   If  DELMAG  is 
zero,  control  is  passed  to  the  garbage  can  lid  section,  since  steepest 
descent  will  not  improve  the  value  of  the  unknown.  If  not,  DELMAG  is  tested 
to  be  certain  it  is  less  than  or  equal  to  one,  and  DELF  is  scaled  down  to 
yield  a  magnitude  of  no  more  than  unity  if  necessary.  This  is  done  to  limit 
the  maximum  change  in  Z  to  a  value  which  will  prevent  crossing  the  "ridge" 
between  two  minima.  A  new  estimate  of  Z  is  obtained  from  an  expression  of 
the  form 

Z  =  ZO  -  DELTA*DELF 
where  ZO  is  the  present  value  of  Z  and  DELTA  is  a  small  number  which  also 
limits  the  maximum  change  in  Z. 

If  this  steepest  descent  step  reduces  the  error,  the  pass  through  the 
iteration  loop  is  complete,  and  the  next  pass  is  begun.   If  the  error  is  not 
reduced,  the  step  distance  DELTA  is  halved,  the  steepest  descent  step 
repeated,  and  the  error  checked  again.  The  step  distance  may  be  halved 
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twenty-six  times  before  steepest  descent  is  abandoned.  This  results  in  a 
minimum  step  distance  which  is  lost  due  to  machine  precision  when  DELTA  and 
Z  are  within  the  usual  limits  of  the  program.  Since  it  is  possible  that 
failure  to  find  a  smaller  error  may  be  due  to  an  incorrect  derivative 
approximation  sign  rather  than  to  having  located  a  minimum,  the  garbage  can 
lid  section  is  entered  after  steepest  descent  fails. 

After  calculating  the  error  at  each  of  twenty-four  points  on  a  circle 
about  the  present  value  of  Z,  the  garbage  can  lid  procedure  obtains  the  min- 
imum error  of  the  set  along  with  the  corresponding  value  of  Z  and  tests  to 
see  whether  or  not  the  error  has  been  improved.  If  the  error  is  reduced, 
the  pass  through  the  iteration  loop  is  complete.   If  not,  the  circle  radius 
is  halved  and  so  on  as  in  the  steepest  descent  portion  of  the  loop.  If  this 
fails  to  find  a  smaller  error,  the  iteration  loop  is  exited. 

Thus,  there  are  three  ways  to  leave  the  iteration  loop:  negligible 
error,  inability  to  find  a  smaller  error,  and  one  hundred  iterations  com- 
pleted without  finding  a  minimum. 

The  IBM  System  360  FORTRAN  IV-G  language  does  not  accept  complex  argu- 
ments in  the  tangent  function,  so  the  function  ZTAN  was  written  to  provide 
this  feature.  Another  problem  arises  when  the  real  part  of  the  argument 
passed  to  ZTAN  is  too  near  an  odd  multiple  of  tt/2;  the  built  in  TAN  function 
will  terminate  the  job.  To  avoid  this,  ZTAN  tests  the  real  part  and  gener- 
ates a  large  number  within  the  range  of  the  machine  to  use  in  place  of  the 
actual  TAN  value  if  the  real  part  of  the  argument  is  near  enough  to  result 
in  job  termination. 

There  are  a  number  of  statements  in  the  basic  program  which  are  not 
described  in  the  preceding  description  or  in  the  flowchart  of  Appendix  A 
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since  they  are  not  essential  to  the  algorithm.  The  variable  IERR  is  used  to 
indicate  whether  or  not  the  argument  proximity  to  tt/2  problem  occurred  in 
ZTAN.  FLAG  is  used  to  trace  the  flow  through  the  loop;  if  FLAG  is  even, 
linear  iteration  did  not  improve  the  error,  if  odd,  linear  iteration  did 
improve  the  error.   If  0  or  1,  steepest  descent  was  used  and  the  garbage  can 
lid  section  was  not;  if  2  or  3,  the  garbage  can  lid  section  was  used  and 
steepest  descent  was  not;  if  4  or  5,  both  were  used.   IFORM  determines  the 
type  of  output  format  used;  this  will  be  described  in  Section  3.3. 

3.3   The  Monitor  and  Control  Section  of  the  Program.  A  flowchart  of 
the  monitor  and  control  section  is  given  in  Appendix  A.  The  program  begins 
with  a  number  of  comment  cards  describing  its  purpose  and  use,  followed  by 
the  non-executable  statements.  After  the  non-executable  statements,  the 
error  function  is  defined  by  an  arithmetic  statement  function.  The  expres- 
sion used  is 

F(Z)  =  Z  -  ER*RATnW0PI*CSQRT((K2SQ-KlSQ)/(ER**2+ZTAN(Z)**2)). 
In  this  expression,  ER  is  the  relative  permittivity  of  region  2  relative  to 
region  1,  RAT  is  the  ratio  of  the  slab  thickness  to  the  free  space  wave- 
length, TWOPI  =  2. *3. 1415927,  and  K2SQ  and  K1SQ  are  symbols  standing  for  the 
product  of  the  relative  permittivity  and  relative  permeability  of  regions  2 
and  1,  respectively.  Following  this  there  are  several  replacement  state- 
ments which  initialize  various  switches  and  variables.  A  page  of  Hollerith 
output  describing  the  output  of  the  program  is  produced  next.  The  entire 
set  of  data  cards  is  read  in  next,  each  card  being  written  on  the  printer  for 
an  echo  check  and  on  a  disk  file  for  re-reading.  When  end-of-file  on  the 
card  read  unit  is  encountered,  the  disk  working  unit  is  rewound;  the  data 
reading  loop  is  then  begun. 
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In  this  data  reading  loop,  a  record  is  read  from  the  disk  working  unit. 
The  format  of  this  READ  statement  calls  for  two  4-character  alphameric 
fields  and  three  10  character  floating  point  numeric  fields.  The  two  alpha- 
meric fields  are  received  by  the  variable  INVARI,  which  is  compared  with 
alphameric  character  strings  stored  in  the  array  LABEL.  When  a  match  is 
found,  the  values  in  the  numeric  fields  are  assigned  to  the  appropriate 
variables.  For  example,  if  the  title  on  the  data  card  is  'FREQ    ',  the 
variable  FGHZ  is  assigned  the  value  of  VALUEA,  the  number  in  the  first  of 
the  three  numeric  fields  on  the  data  card.  After  such  a  match  is  found  and 
the  appropriate  value  or  values  assigned  to  the  appropriate  variable  or 
variables,  another  card  is  read  from  the  disk  working  unit.  The  last  card 
of  the  deck  (the  last  FORTRAN  data  card)  should  be  a  card  containing  the 
character  string  'STOP    ';  this  card  causes  subroutine  EXIT  to  be  executed, 
If  no  match  is  found,  the  erroneous  character  string  is  printed  out,  fol- 
lowed by  a  description  of  the  allowable  card  types  and  formats;  then  execu- 
tion is  terminated  by  the  statement  STOP. 

The  appearance  of  a  card  'RUN  '  causes  the  remainder  of  the  program 
to  be  executed  using  the  parameters  in  effect  at  that  time;  after  such  exe- 
cution, another  card  is  read. 

Thus,  the  program  may  produce  several  sets  of  output  with  a  minimum  of 
input;  if  all  that  is  to  be  changed  between  two  sets  of  parameters  during 
the  same  job  is  the  relative  permittivity  of  the  dielectric  (region  2),  the 
data  cards  following  the  first  run  card  (that  is,  the  data  cards  for  the 
second  set  of  parameters)  need  only  be  1 )  a  card  containing  the  title 
' ER2P    '  and  the  new  value  in  the  first  numeric  field,  and  2)  the  second 
'RUN    '  card  (followed  by  more  data  sets  or  the  'STOP     '  card).  For  a 
reasonably  detailed  description  of  the  allowable  data  cards,  refer  to  the 
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comment  cards  in  the  program  listing  which  is  contained  in  Appendix  B.  Some 
of  the  cards,  which  serve  only  control  functions  such  as  do  the  'RUN     ' 
and  'STOP     '  cards,  contain  only  a  title;  others  contain  one  or  two  nu- 
meric values.  A  notable  exception  is  $DATA,  the  array  used  to  set  the  value 
of  the  loss  component  of  permittivity  or  permeability  during  loss  parameter 
sweeps:  the  cards  to  determine  the  values  of  the  elements  of  $DATA  and  the 
number  of  active  elements  are  1)  a  card  containing  the  title  '$DATA    '  and 
the  number  of  active  elements  in  the  first  numeric  field,  and  2)  cards  con- 
taining the  values  of  the  elements  in  10F8.4  format.  For  the  mode  numbers, 
the  integers  should  be  even  since  the  program  is  valid  only  for  even  mode 
numbers.  Note  that  all  the  numeric  values  on  the  data  cards  are  floating 
point,  and  that  some  are  converted  to  integer  values.  The  statement  "GO  TO 
900"  preceding  statement  924  is  the  last  statement  of  the  data  reading  loop. 

The  appearance  of  a  'RUN     '  card  causes  control  to  be  passed  to 
statement  924.  The  parameters  ER,  K1SQ,  and  K2SQ  are  calculated.  Next  the 
mode  loop  is  entered;  this  is  simply  a  DO  loop  having  index  parameters 
related  to  the  lowest  and  highest  mode  numbers  requested.  The  variable  IM 
contains  the  present  mode  number  throughout  the  mode  loop. 

Next  the  starting  value  of  Z  is  determined.   If  a  starting  value  has 
been  given  as  input  data,  as  indicated  by  the  variable  IZSW  having  a  value 
of  1 ,  Z  is  set  equal  to  the  given  value;  control  is  passed  to  statement  925, 
bypassing  the  remainder  of  the  starting  value  calculations.   If  IZSW  is 
zero,  the  program  calculates  a  starting  value  for  Z  by  the  following  proce- 
dure: the  radius  of  the  circle  (with  reference  to  Figure  2)  is  calculated 
using  the  real  parts  of  the  relative  permittivities  and  permeabilities.  If 
the  radius  is  greater  than  or  equal  to  IM*PI/2,  the  mode  is  assumed  to  be  a 
surface  wave  (propagating)  mode  and  control  is  passed  to  statement  955.   If 
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this  radius  condition  is  not  met,  the  mode  is  assumed  to  be  a  leaky  wave 
(cutoff)  mode  and  Z  is  calculated  according  to  the  statement 

Z  =  (IM  -  1  +  0.1)*PI/2  +  J*0.5 
where  J  stands  for  the  imaginary  unit  and  PI  =  3.1415927;  control  is  passed 
to  statement  925.  At  statement  955,  the  radius  is  tested  to  determine 
whether  or  not  it  is  greater  than  or  equal  to  (IM+l)*PI/2.   If  not,  then  it 
must  be  between  IM*PI/2  and  (IM+l)*PI/2,  so  a  starting  value  in  this  range 
is  determined  from  the  statement 

Z  =  0.9*RADIUS  +  0.1*IM*PI/2; 
control  is  transferred  to  statement  925.  If  the  test  expression  is  true, 
control  is  transferred  to  statement  960  where  Z  is  calculated  from  the 
expression 

Z  =  (IM+0.9)*PI/2. 
Statement  925  follows  statement  960. 

NPASS  is  a  variable  indicating  the  number  of  times  the  iteration  loop 
has  been  executed.  IELEC  and  IMAG  are,  respectively,  indices  which  deter- 
mine which  element  of  $DATA  is  used  during  electric  and  magnetic  sweeps  for 
ER2PP  or  MUR2PP,  the  loss  components  of  the  relative  permittivity  and  rela- 
tive permeability  of  region  2.   IESW  determines  whether  or  not  an  electric 
sweep  is  requested,  while  IMSW  performs  the  same  function  for  magnetic 
sweeps.   If  both  sweeps  are  requested,  the  electric  sweep  will  be  performed 
first.  After  the  appropriate  parameters  have  been  modified,  the  iteration 
loop  is  entered.  The  value  of  Z  having  thus  been  determined,  U,  V,  KZ, 
LRAT,  ATTX,  and  ATTZ  are  calculated.  U  is  the  x-directed  component  of  the 
propagation  constant  in  region  2,  V  is  the  x-directed  component  of  the 
propagation  in  region  1,  KZ  is  the  z-directed  component  of  the  propagation 
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constant  in  both  regions,  ATTX  is  the  attenuation  in  the  x-direction  in 
region  1,  and  ATTZ  is  the  attenuation  in  the  z-direction;  all  these  param- 
eters are  in  nepers  per  free-space  wavelength  except  ATTX  and  ATTZ,  which 
are  in  decibels  per  free-space  wavelength.  LRAT  is  the  ratio  of  the  free- 
space  wavelength  to  the  guide  wavelength.  The  output  is  written  next.  The 
format  of  the  output  is  determined  by  the  variable  IFORM.  IFORM  =  0  pro- 
duces an  output  tracing  the  value  of  Z,  the  error,  and  the  derivatives 
through  the  iterations  and  allows  messages  warning  of  the  use  of  the  approx- 
imation within  ZTAN  to  be  printed.  This  output  is  suitable  for  checking  the 
proper  operation  of  the  program.   IFORM  =  1  produces  a  tabular  output  of  the 
results  only,  giving  no  indication  of  the  magnitude  of  the  error,  the  number 
of  iterations,  or  other  such  information.  After  producing  the  output  for  a 
pass,  NPASS  is  checked;  on  the  first  pass,  Z  is  stored  in  ZSTART  for  use  in 
the  magnetic  sweep,  should  one  be  requested.  IELEC  and  IMAG  are  checked 
next;  if  one  of  the  sweeps  is  finished,  the  appropriate  parameter  ER2PP  or 
MUR2PP  is  reset  to  zero.  If  an  electric  sweep  is  in  progress  but  incomplete, 
IELEC  is  incremented  and  the  next  pass  performed;  similarly  for  a  magnetic 
sweep  IMAG  is  incremented.  After  completing  the  sweeps  or  the  single  pass 
if  no  sweeps  are  requested,  the  process  is  repeated  for  the  next  mode.  When 
the  mode  loop  is  satisfied,  control  is  passed  to  the  READ  statement  number 
900  to  begin  the  input  reading  loop  again. 

3.4   Comments  on  Programming.  The  emphasis  in  this  program  has  been 
placed  strictly  upon  convenience  to  the  programmer,  with  no  attention  paid 
to  the  execution  time.  For  this  reason,  the  symbol  J  was  defined  to  be  the 
complex  constant  (0.,  1.);  complex  expressions  were  then  written  in  the  form 
A  +  J*B.  This  requires  the  conversion  of  A  and  B  to  complex  mode,  a  complex 
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multiplication,  and  a  complex  addition.  The  same  result  can  be  obtained  by 
the  use  of  the  in-line  built-in  function  CMPLX(A,B)  with  a  very  considerable 
reduction  in  execution  time.  Should  the  program  be  used  for  further  re- 
search into  a  large  number  of  cases,  it  would  be  highly  appropriate  to  re- 
place the  statements  using  J  with  statements  referencing  CMPLX(A,B). 

There  are  several  peculiarities  in  the  program  which  deserve  discussion. 
The  most  important  of  these  is  the  starting  value  used  for  Z  and  its  rela- 
tion to  the  mode  number  to  which  the  result  corresponds.  For  the  lossless 
surface  wave  case,  the  value  of  Z  is  real  and  clearly  bounded  between 
IM*PI/2  and  (IM+l)*PI/2.  Thus,  a  starting  value  between  these  bounds  is 
virtually  certain  to  yield  the  desired  result.  For  the  other  cases,  the 
root  cannot  easily  be  bounded,  so  that  a  simple  expression  for  a  suitable 
starting  value  has  not  been  determined.  It  is  quite  possible  that,  from  a 
given  starting  value,  the  algorithm  may  converge  to  a  root  corresponding  to 
a  different  mode  than  that  which  the  programmer  specified.  It  is  necessary 
for  the  programmer  to  determine  whether  or  not  the  root  does  correspond, 
and  if  necessary,  to  insert  an  appropriate  starting  value  by  means  of  a  data 
card. 

The  Cutoff  modes  are  particularly  troublesome  in  this  respect.  Since, 
in  the  lossless  case,  the  error  is  an  even  function  of  the  imaginary  part  of 
Z,  a  real  starting  value  must  not  be  used,  as  the  algorithm  in  this  event 
will  converge  to  a  non-zero  minimum  on  the  real  Z  axis.  For  severely  cutoff 
modes  (RADIUS  <<  IM*PI/2),  the  root  may  be  so  near  an  odd  multiple  of  tt/2 
that  it  is  necessary  to  specify  a  small  value  of  DELTAO  to  prevent  the  algo- 
rithm from  allowing  Z  to  jump  over  a  "wall"  in  the  error  space. 

When  the  sweep  feature  is  used  and  the  first  case  is  a  lossless  surface 
wave  case,  the  mode  uncertainty  problem  is  not  likely  to  occur.  That  is, 
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the  expressions  used  to  calculate  a  starting  value  of  Z  generate  a  value  in 
the  proper  interval.  The  program  keeps  the  last  root  value  as  the  starting 
value  for  the  following  case  during  a  sweep,  so  appealing  to  some  notion  of 
continuity,  the  correct  roots  should  be  obtained  as  the  losses  are  increased 
if  the  increments  of  the  swept  variable  are  not  too  large.  For  lossless 
cutoff  cases  in  which  the  degree  of  cutoff  is  not  too  severe,  the  built-in 
determination  of  the  starting  value  should  be  satisfactory. 

The  garbage  can  lid  section  was  included  in  order  to  cope  with  the 
problem  of  a  calculated  gradient  of  zero  when  in  fact  the  value  of  Z  was  not 
sufficiently  close  to  a  zero  of  the  error.  It  was  necessary  to  iterate  un- 
til the  root  was  accurate  to  7+  significant  digits  in  order  to  obtain  suit- 
ably accurate  values  of  KZ  and  V.  This  may  in  part  be  due  to  the  calcula- 
tion of  KZ  from  the  expression 

KZ  =  CSQRT(4.*PI*PI*(MUR2P-J*MUR2PP)*(ER2P-J*ER2PP)-U**2) 
which  results  in  taking  the  difference  of  nearly  equal  imaginary  parts  for 
lossy  cases.  Later  versions  of  the  program  used  the  expression 

KZ  =  CSQRT(4.*PI*PI*(MUR1P-J*MUR1PP)*(ER1P-J*ER1PP)+V**2) 
which  gives  improved  accuracy  particularly  when  region  1  is  lossless.  It  is 
still  necessary  to  obtain  the  maximum  single  precision  accuracy  in  the  root 
for  large  relative  permittivity  cases. 


CHAPTER  IV 
INTERPRETATION  OF  DATA 

4.1  Behavior  of  TMo  Surface  Wave  Modes.  The  TMo  mode  always  propa- 
gates as  a  surface  wave  and  is  the  dominant  mode.  It  is  probably  the  most 
interesting  mode  for  practical  purposes. 

The  results  from  the  computer  analysis  for  the  TMo  mode  are  presented 
in  graphical  form  in  Figures  5  through  18;  reproductions  of  the  computer 
output  are  contained  in  Appendix  C.  Figures  5  and  6  show  the  z-direction 
attenuation  as  a  function  of  e  «"  for  e  «'  =  2  and  values  of  (t/Xo)  from 
.005  to  .50. 

The  same  results  may  be  obtained  in  closed  form  for  t/Xo  <  .05  or  so  by 
approximating  tan  (ut)  by  (ut)  in  Equation  2.1.18  and  solving  the  resulting 
quartic  for  u.  This  has  been  done  for  the  low  loss  case  by  Barlow  and 
Cullen  [1953]  and  Johnson  [1966].  As  might  be  expected,  this  approximation 
gives  a  significant  error  at  t/Xo  =  .10  and  does  not  predict  the  interesting 
results  which  occur  between  t/Xo  =  .10  and  t/Xo  =  .20. 

In  general,  the  attenuation  tends  to  increase  as  t/Xo  increases  up  to 
t/Xo  =  .10.  For  t/Xo  =  .18  and  below  the  curves  exhibit  a  tendency  to  turn 
downward  at  high  losses  while  for  t/Xo  =  .19  and  above  the  curves  exhibit  a 
tendency  to  turn  upward.   In  this  transition  region  the  curves  break  sharply, 
with  the  actual  switch  occurring  between  t/Xo  =  .18  and  t/Xo  =  .19,  and  at  an 
e  2     of  approximately  2.  As  the  ratio  increases  from  .10  to  .75  the  curves 
indicate  that  for  small  e  «"  the  attenuation  increases  as  t/Xo  increases  but 
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Figure  5.  z-di recti  on  Attenuation  for  the  TMo  Mode  and 
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Figure  6.  z-direction  Attenuation  for  the  TMo  Mode  and  e  ~'=2  (Continued) 
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for  large  e  -"  the  attenuation  decreases  as  t/Xo  increases.  That  is,  com- 
paring two  curves,  the  curve  for  the  larger  ratio  will  be  above  the  curve 
for  the  smaller  ratio  when  e  «"  is  small  but  the  curves  will  cross  at  s 
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value  of  e  «"  (in  the  vicinity  of  e."  =  2  f°r  t/Xo  >  .19)  and  reverse  rela- 
tive positions  for  larger  e  «"  values.  Obviously  this  is  not  true  for  a 
comparison  of  ratios  of  .18  and  .19,  but  it  is  true  if  the  two  values  under 
consideration  are  both  between  .10  and  .18  or  both  greater  than  .19.  For 
ratios  greater  than  .50,  at  which  value  the  second  mode  becomes  feasible  as 
a  surface  wave,  there  is  little  change  in  the  z-di recti  on  attenuation  as 
t/Xo  is  increased. 

Figure  7  shows  the  x-direction  attenuation  for  the  same  cases  as  in 
Figure  5.  There  is  a  great  deal  of  similarity  between  these  curves  and  the 
curves  of  z-direction  attenuation;  they  are  related  through  Equation  2.1.15, 
the  separation  equation  in  region  1.  Since  the  attenuation  involves  the 
real  part  of  v  and  the  imaginary  part  of  k  ,  the  relation  is  not  easy  to 
visualize  from  the  formula.  Some  curves  representative  of  the  guide  wave- 
length along  the  z  axis  are  shown  in  Figure  8.  In  the  lossless  case,  the 
guide  wavelength  is  always  shorter  than  the  free  space  wavelength  and  hence 
the  phase  velocity  is  less  than  the  speed  of  light  in  free  space.  As  losses 
are  inserted,  the  guide  wavelength  decreases  slightly  for  yery   small  t/Xo 
but  for  t/Xo  values  in  the  transition  range  initially  increases  in  a  yery 
pronounced  fashion,  peaks  for  e  «"  near  2  as  do  the  attenuation  curves,  and 
drops  off  for  further  increase  in  e  p"-  In  some  cases  the  guide  wavelength 
and  phase  velocity  rise  above  their  free-space  values. 

Figures  9,  10,  and  11  for  e  «'  =  4.24  correspond  to  Figures  5,  7,  and 
8  respectively;  similarly  the  case  for  e  2'  =  8.28  is  presented  in  Figures 
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Figure  7. 


x-direction  Attenuation  for  the  TMo  Mode  and  e  ?'=2 
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Figure  9.  z-direction  Attenuation  for  the  TMo  Mode  and  e  0 ' =4 .24. 
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Figure  10.  x-di recti  on  Attenuation  for  the  TMo  Mode  and  e  ?'=4.24 
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12,  13,  and  14,  while  the  case  for  e  2'  =  165  is  presented  in  Figures  15, 
16,  17,  and  18.  The  features  of  these  families  of  curves  are  essentially 
identical  to  those  for  e  ?'  =  2.  Table  2  shows  the  approximate  values  of 
e  p"  and  t/Ao  at  which  the  transition  occurs  for  each  value  of  e  ?'. 

TABLE  2 
APPROXIMATE  er2"  AND  t/Ao  TRANSITION  VALUES 


er2' 

2 

4.24 

8.28 

165 

11 
er2 

2 

3 

4.1 

18 

t/Ao 

.18 

.115 

.083 

.018 

Roughly  speaking,  the  following  empirical  relations  among  the  transition 
values  hold: 


er2M  =  /^7  [4'1J] 


t/Ao  *  .24/  /~/    .  [4.1.2] 


Mathematically,  it  appears  that  the  transition  phenomenon  is  due  to  behavior 
of  tan  (z)  for  Re(z)  near  tt/2;  at  the  transition,  the  value  of  Re(z)  is  be- 
tween .43tt  and  .48tt.  Some  further  physical  and  mathematical  insight  as  to 
the  mechanism  of  transition  will  be  gained  in  later  paragraphs. 

In  section  2.3  it  was  noted  that  conjugate  pair  solutions  occur  for  the 
lossless  leaky  wave  case,  and  it  was  conjectured  that  approximately  conju- 
gate roots  might  exist  in  the  lossy  surface  wave  case.  Such  roots  have  been 
encountered.  For  cases  in  which  region  1  is  lossy,  such  roots  may  have 
physical  significance;  they  have  not  been  examined  in  this  research. 
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Figure  12.  z-di recti  on  Attenuation  for  the  TMo  Mode  and  e  ?'=8 


=8.28. 
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x-direction  Attenuation  for  the  TMo  Mode  and  e  0,=8.28. 
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Figure  14.  Wavelength  Ratio  for  the  TMo  Mode  and  e  2'=8.28. 
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Figure  15.  z-direction  Attenuation  for  the  TMo  Mode  and  e  2 ' =^ 65 . 
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Figure  16.  z-di recti  on  Attenuation  for  the  TMo  Mode  and  e  2'=165  (Continued) 
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Figure  17.  x-di  recti  on  Attenuation  for  the  TMo  Mode  and  e^'3^. 
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Figure  18.  Wavelength  Ratio  for  the  TMo  Mode  and  e  ?'=165, 


50 


Inspection  of  the  computer  output  shows  that,  as  losses  are  inserted, 

e„  enters  the  fourth  quadrant  while  e  enters  the  first  quadrant.  That  is, 

P  a 

the  direction  of  propagation  tilts  so  as  to  carry  energy  into  the  slab  from 
without;  one  may  hypothesize  that  this  tilt  occurs  in  order  to  overcome  the 
losses  of  the  dielectric.  If  region  1  is  lossless,  as  is  true  in  all  cases 
considered  in  this  research,  the  wave  in  region  1  is  still  an  evanescent 

wave,  that  is  ore  =  0  or  e  =  8  +  tt/2 ;  this  is  an  immediate  consequence  of 

a    p 

Equation  2.1.15.  The  angle  of  incidence  is  tt/2  +  e  =  e  .  As  might  be 

expected  from  the  data  already  presented,  e.  is  zero  in  the  lossless  case, 

p 

becoming  more  and  more  negative  with  the  insertion  of  losses  until  the  tran- 
sition effect  sets  in  (see  Figure  21).  For  t/Xo  less  than  the  transition 
value,  e  reaches  a  maximum  negative  excursion  for  some  value  of  e  «"  and 
then  tends  again  to  parallel  incidence  as  losses  are  further  increased.  For 

t/Ao  greater  than  the  transition  value  eD  appears  to  tend  toward  a  constant 

p 

which  tends  to  -tt/4  as  t/Xo  increases.  This  may  be  the  best  explanation  of 
decreasing  attenuation  with  increasing  e  a"* 

Though  a  small  amount  of  research  into  the  behavior  of  dielectrics 
having  magnetic  losses  has  been  performed,  sufficient  data  has  not  been 
collected  to  discuss  the  behavior  of  such  cases.  Their  behavior  can  be 
expected  to  be  slightly  different  since  the  e  term  is  not  affected  by  mag- 
netic losses. 

4.2   The  Behavior  of  the  TM-.  Extraordinary  Surface  Wave  Mode.  When 
the  transition  effect  was  first  encountered,  it  was  suspected  that  there  was 
some  mathematical  error  in  the  program,  such  as  an  improper  choice  of  sign 
preceding  a  square  root.  A  hand  check  indicated,  however,  that  the  values 
obtained  were  indeed  valid  roots  of  Equation  2.1.18.  The  value  of  u  (per 
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free  space  wavelength)  was  plotted  for  the  various  values  of  t/Xo;  e  ~'  =   2 
was  being  investigated  at  the  time  in  question.  This  plot  is  included  as 
Figure  19.  The  curves  for  the  TMo  ordinary  mode,  i.e.,  those  originating 
on  the  Real  u  axis,  were  the  only  ones  plotted  at  this  point  in  time.  A  run 
was  made  for  t/Xo  =  .19  and  .20  for  an  initial  e  2"  of  2.2  using  a  ZSTART 
near  the  U*RAT  value  previously  produced  by  t/Xo  =  .18  with  the  intention  of 
verifying  that  the  value  of  u  would  indeed  meander  away  from  this  starting 
value  to  the  values  previously  obtained  for  t/Xo  =  .19  and  .20.  The  DEBUG 
format  option  was  specified  to  obtain  visibility  into  the  path  of  conver- 
gence and  the  error  magnitude.  The  roots  which  were  found,  however,  were 
not  the  same  as  before--they  were  quite  near  the  u  values  for  t/Xo  =  .18. 
What  a  revolting  development.  The  next  logical  step  was  to  sweep  e  «"  d°wn 
to  0  and  up  to  6  and  plot  these  results  in  order  to  obtain  further  insight. 
This  was  done  for  ratios  of  .19  and  .20;  also,  for  t/Xo  =  .18  and  the  single 
point  e  «"  =  6>  t/Xo  =  .16  solutions  were  anticipated  corresponding  to  the 
ordinary  mode  behavior  for  t/Xo  =  .19  and  .20.  All  these  values  were  indeed 
found;  indicative  curves  of  u  values  for  this  extraordinary  mode  are  shown 
in  Figure  19. 

The  curves  extending  below  Re(u)  =  8.75  are  those  of  the  ordinary  mode 
while  those  not  extending  below  Re(u)  =  8.75  correspond  to  the  extraordinary 
mode.  The  boundary  separating  the  ordinary  and  extraordinary  modes  has 
been  empirically  determined.  Those  having  decreasing  Im(u)  at  high  e  2" 
produce  decreasing  z-direction  attenuation  at  high  e  J1   while  those  having 
increasing  Im(u)  at  high  e  2"  produce  increasing  z-direction  attenuation  at 
high  er2". 

Since  the  DEBUG  format  option  was  used  during  this  investigation,  the 
computer  output  is  too  bulky  to  include.  However,  tabular  output  for  the 
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extraordinary  mode  with  e  «'  =  165>  t/Xo  =  .02  is  included  in  Appendix  C 
(the  mode  subscript  is  incorrect,  if  indeed  it  is  meaningful). 

It  is  now  apparent  that  this  extraordinary  mode  corresponds  to  a  mode 
numbered  TM-, .  Such  modes  are  not  physically  real  in  the  lossless  case  where 
there  is  a  real  intersection  of  the  z-tan(z)  curve  with  the  circle  of  radius 


2-rr[t/Ao]/y  ~z   2  -   y  -, e  , ;  3  lies  in  the  a  direction  but  a  lies  in  the  -a 
direction.  Examination  of  the  computer  results  shows  that  this  TM-,  mode  has 
an  a  lying  in  the  third  quadrant  or  on  its  boundary  for  low  or  zero  losses, 
but  as  losses  are  increased  a  passes  through  zero  magnitude  and  enters  the 
first  quadrant;  0  remains  in  the  fourth  quadrant  in  both  cases.  The  zero 
magnitude  of  a  is  easy  to  believe  for  the  case  shown  in  Appendix  C  but  may 
be  an  erroneous  conclusion  for  the  general  case.  That  is,  from  the  require- 
ment k  * k  "  =  v'v",  obtained  from  Equation  2.1.15  with  k,  real,  and  the 
smooth  transition  of  v'  through  zero  it  is  apparent  that  k  '  k  "  passes 
through  zero,  k  "  changes  sign;  k  '  does  not  change  sign  but  may  approach 
zero,  possibly  allowing  k  "  to  be  discontinuous.  A  closer  investigation  of 
the  region  1.4  <  e  2"  <   I-6  should  be  made  for  e  2'  =  2,  t/Xo  =  .18  and  .19. 

The  absolute  value  of  the  z-direction  attenuation  for  the  TM-,  extra- 
ordinary mode  is  shown  in  Figure  20.  The  dotted  sections  are  actually  nega- 
tive, the  change  in  direction  of  a  occurring  between  e  ?"  =  1.4  and 
e  2     =   1'6»  The  transition  effect  also  appears  in  this  mode.  In  fact,  the 
extraordinary  mode  for  t/Xo  >_  .19  continues  the  trend  of  the  TMo  surface 
wave  mode  for  t/Xo  ^  .18  and  the  extraordinary  mode  for  t/Ao  <_  .18  continues 
the  trend  of  the  TMo  wave  mode  for  t/Xo  >_  .19.  Figure  21  shows  this  comple- 
mentary relationship  for  the  angle  e  . 

4.3   Behavior  of  the  TMq  Leaky  Wave  Mode.  For  a  circle  having  a  radius 
less  than  tt,  the  TM2  mode  is  a  cutoff  mode  giving  rise  to  a  leaky  wave,  t 
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Figure  20.     z-direction  Attenuation  for  the  TM,  Mode  and  erp,=2 
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lies  in  the  first  quadrant  and  a  in  the  fourth  quadrant  when  region  1  is 
lossless. 

Some  representative  computer  output  is  included  in  Appendix  C.  Only 
cases  in  which  region  2  is  a  dielectric  having  e  ?'  =  2  with  the  magnetic 
properties  of  free  space  and  region  1  has  the  properties  of  free  space  were 
investigated.  The  z-direction  attenuation  is  not  of  as  much  interest  here 
as  in  the  surface  wave  cases,  so  it  is  not  plotted.  In  general,  z-direction 
attenuation  increases  with  decreasing  t/Xo,  i.e.,  a  more  severe  degree  of 
cutoff,  e  is  a  much  more  interesting  quantity.  This  angle,  which  might 
appropriately  be  termed  the  angle  of  departure  or  launch  angle  in  the  leaky 

wave  cases,  is  plotted  in  Figure  22.  e0  increases  as  the  degree  of  cutoff 

p 

increases.  For  t/Xo  =  .49  the  z-direction  attenuation  and  launch  angle  are 
very   nearly  zero  for  e  «"  =  0  since  the  degree  of  cutoff  is  very  mild. 
There  may  be  effects  of  interest  occurring  in  the  range  .10  <  t/Xo  <  .30 
which  are  not  shown  as  evidenced  by  the  dissimilarity  of  the  t/Xo  =  .20 
curve  to  the  t/Xo  =  .10  and  .30  curves.  The  wave  seems  to  be  affected  very 
little  by  losses  in  severely  cutoff  cases. 

A  few  notes  are  necessary  concerning  the  cases  t/Xo  =  .49  and 
t/Xo  =  .01.  For  e  2"  =  0  and  t/Xo  =  .49  the  computer  oscillated  between  the 
proper  solution  and  the  conjugate  solution,  finally  converging  to  the  con- 
jugate solution.  This  explains  the  sign  discrepancy;  format  field  width 
explains  the  0.00000  value:  the  imaginary  part  of  z  is  only  4.5E-10,  thus 
the  two  solutions  are  extremely  close  together.  It  is  simply  fortuitous  that 
the  error  surface  is  usually  inclined  such  as  to  lead  the  iteration  routine 
to  the  proper  solution;  no  planning  or  constraints  forced  it.  Only  the  first 
line  of  output  is  valid  for  t/Xo  =  .01;  the  degree  of  cutoff  is  extreme  so 
that  Re(z)  is  near  tt/2:   1.5710928  versus  1.5707963.  This  is  an  example  of 
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the  increment  magnitude  problem  mentioned  in  Section  3.4.  At  the  time  of 
this  writing,  attempts  to  obtain  valid  results  for  t/Xo  =  .01  have  been 
unsuccessful . 

Some  conjecture  about  the  possibility  of  an  extraordinary  cutoff  mode 
seems  appropriate.  The  conjugate  solution  to  the  TM-,  extraordinary  mode 
could  exhibit  the  characteristics  of  a  leaky  wave;  referring  to  the  TM-, 
extraordinary  mode  e  ?'  =  165  computer  output  in  Appendix  C,  a  conjugate 
solution  would,  at  low  losses,  give  rise  to  a  I  in  the  first  quadrant  and 
an  a  in  the  fourth  quadrant. 

4.4   Behavior  of  the  TMq  Surface  Wave  Mode.  The  case  of  e  «'  =  2  was 
investigated  for  the  TNL  surface  wave  mode,  the  properties  of  region  1  and 
the  magnetic  properties  of  region  2  being  those  of  free  space.  The  z-direc- 
tion  attenuation,  x-direction  attenuation,  and  wavelength  ratio  are  plotted 
in  Figures  23,  24,  and  25  respectively.  The  behavior  of  this  mode  is  highly 
similar  to  that  of  the  TMo  mode.  There  are  three  notable  differences: 
first,  the  value  of  t/Ao  at  which  the  transition  occurs  is  different;  second, 
the  value  of  e  ~"   at  which  the  transition  occurs  is  different;  third,  for 
high  losses  and  t/Xo  less  than  the  transition  value  the  curves  for  various 
t/Xo  values  are  so  close  together  that  it  is  impossible  to  distinguish  one 
from  another. 

One  might  anticipate  a  difference  in  the  TMo  and  TM~  modes  for  circle 
radii  near  the  cutoff  value  in  each  case  since  the  z-tan(z)  curve  has  zero 
slope  at  z  =  0  but  non-zero  slope  at  z  =  tt.  Such  a  difference  is  not  appar- 
ent, however.  On  the  other  hand,  compare  the  region  0  <  e  «"  <  2  for  the 
TM?  with  the  region  0  <  e  2"  <  6  for  the  TMo;  the  TMo  might  exhibit  the  same 
asymptotic  behavior  as  the  TMp  if  e  «"  were  increased  beyond  6.  This  might 
be  the  manifestation  of  the  slope  difference. 
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Figure  23 


z-di recti  on  Attenuation  for  the  TMp  Surface  Wave  Mode  and  e ro'=2 
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Figure  24, 


x-di recti  on  Attenuation  for  the  TNL  Surface  Wave  Mode  and  e  0'=2 
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Figure  25.  Wavelength  Ratio  for  the  TMp  Surface  Wave  Mode  and  e  ? 
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CHAPTER  V 
SUMMARY  AND  CONCLUSIONS 

5.1   Summary.  Chapter  II  begins  with  the  derivation  of  the  charac- 
teristic equation  and  the  elimination  of  variables  to  obtain  a  single  equa- 
tion in  one  unknown.  A  graphical  analysis  of  the  real  roots  is  presented; 
this  graphical  solution  has  been  presented  by  several  other  authors  [Collin, 
1960;  Harrington,  1961]  and  is  not  a  significant  contribution.  Chapter  II 
is  concluded  with  a  discussion  of  the  location  and  nature  of  complex  zeros 
of  Equation  2.1.18. 

Chapter  III  begins  with  a  discussion  of  the  numerical  approaches  with 
which  early  experiments  were  carried  out.  The  final  scheme  of  iteration,  a 
combination  of  steepest  descent  and  linear  iteration,  is  described  and  the 
FORTRAN  IV  implementation  explained.  The  remainder  of  the  program  is  de- 
scribed in  order  that  other  researchers  may  understand  the  details  of  use  of 
the  program.  The  chapter  is  concluded  with  a  discussion  of  problems  and 
limitations  within  the  program  and  recommendations  for  improving  upon  them. 

Chapter  IV  presents  the  data  compiled  and  its  interpretation.  The 
quantities  of  interest  are  the  attenuation  in  directions  parallel  and  normal 
to  the  slab  (the  z  and  x  directions,  respectively),  the  guide  wavelength,  and 
the  tilt  of  the  wave  outside  the  slab.  These  quantities  are  functions  of  the 
properties  of  both  regions  and  the  thickness  of  the  slab  in  relation  to  the 
free-space  wavelength.  Though  the  computer  program  was  written  to  provide 
for  more  general  properties,  the  electric  and  magnetic  properties  of  the 
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region  above  the  slab  and  the  magnetic  properties  of  the  slab  were  chosen 
to  be  the  properties  of  free  space  in  all  cases  presented  in  this  research. 
The  TMo  mode,  which  is  the  dominant  mode,  shows  the  expected  decrease 
in  z-direction  attenuation  at  high  losses  for  small  thicknesses.  Another 
more  interesting  effect  occurs  as  the  thickness  is  increased;  this  is  the 
sudden  and  drastic  transition  from  decreasing  to  increasing  attenuation  at 
high  losses.  Comparison  of  families  of  curves  for  different  e  «'  values 
indicates  that  for  small  t/Xo,  the  higher  the  e  2'  the  lower  the  attenuation 
while  for  large  t/Xo,  the  higher  the  e   ~'  the  higher  the  attenuation.  The 
case  e  yl   =  1 65  does  not  fit  this  pattern;  the  transition  effects  sets  in  at 
a  much  lower  value  of  t/Ao  than  in  the  other  cases.  Empirically,  the  rela- 
tionship between  the  parameters  at  which  transition  occurs  is  given  by  the 
expressions 


er2     ^er2 
t/Xo  -  .24/  /e 


er2 


Except  for  t/Xo  values  in  the  transition  range,  the  attenuation  in  the  x 
direction  tends  to  increase  slightly  as  losses  are  inserted  and,  as  in  the 
lossless  case,  tends  to  increase  as  t/Ao  increases.  That  is,  as  losses  are 
added,  a  increases  in  magnitude,  since  both  components  increase.  3  tends  to 
tilt  into  the  slab  as  losses  are  added,  the  degree  of  tilt  increasing  with 
increasing  losses  except  near  the  transition  t/Ao  values. 

The  behavior  of  the  wave  for  t/Ao  in  the  transition  range  changes 
drastically.  For  t/Ao  slightly  less  than  the  transition  value,  the  wave 
first  tilts  into  the  slab  as  losses  are  inserted,  then  tilts  back  toward 
parallel  incidence  at  high  losses,  a  increases  at  low  losses  but  decreases 
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at  high  losses  also.  The  picture  is  that  of  a  sufficiently  high  loss  slab 
reflecting  most  of  the  wave  so  that  it  is  less  tightly  bound  and  less  energy 
is  dissipated  in  the  slab.  For  t/Ao  greater  than  the  transition  value,  the 
wave  tilts  into  the  slab,  reaches  some  limiting  angle,  and  continues  to  plow 
into  the  slab  at  a  nearly  constant  angle  as  losses  increase  further;  a  tends 
to  increase  in  magnitude.  Here  the  picture  is  that  the  thicker  slab  prevents 
the  angle  of  incidence  from  returning  toward  zero  so  that  a  great  deal  of 
energy  is  concentrated  in  the  slab  and  dissipated  making  the  wave  tightly 
bound. 

It  was  discovered  that,  for  sufficiently  high  losses,  the  TM, -numbered 
solution  has  the  proper  attributes  to  be  a  surface  wave.  The  behavior  of 
this  extraordinary  mode  is  concisely  but  completely  described  by  saying  that 
it  is  complementary  to  the  TMo  mode. 

Suppose  that  a  source  is  designed  to  excite  surface  waves  on  some  struc- 
ture, and  that  the  value  of  t/Xo  is  larger  than  the  transition  value  but 
small  enough  that  the  TM?  and  higher  modes  are  cut  off.  The  load  impedance 
which  the  source  sees  depends  upon  the  £  and  H  of  each  mode.  The  z  directed 
wave  impedance  tends  to  be  more  and  more  reactive  as  e,  y,  v,  and  k  become 
more  nearly  pure  imaginary.  One  would  expect,  then,  that  the  TM,  mode  would 
carry  most  of  the  power  in  the  structure,  and  that  the  TMo  would  be  difficult 
to  detect.  As  frequency  is  increased  from  a  small  t/Xo  through  transition, 
one  might  expect  a  reasonably  smooth  transition  from  the  TMo  mode  to  the  TM, 
mode  rather  than  an  abrupt  jump  in  attenuation. 

The  TM?  leaky  wave  mode  does  not  appear  to  exhibit  any  particularly 
noteworthy  characteristics.   It  is  expected  that  a  solution  complementary  to 
the  TMp  leaky  wave  mode  may  arise  from  the  extraordinary  TM,  mode. 
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The  TM?  surface  wave  mode  exhibits  properties  very   similar  to  those  of 
the  TMo  mode.  The  merging  of  the  curves  at  high  losses  might  also  occur  for 
the  TMo  mode  at  higher  losses. 

5.2   Recommendations  for  Further  Study.  The  geometry  shown  in  Fig- 
ure 1  will  also  support  TE  modes  where  n  must  be  odd  in  the  lossless  case. 
The  geometry  having  a  mirror  image  half-slab  and  surrounding  medium  in  place 
of  the  grounded  conductor  will  support  TE  and  TNI  modes  having  both  even  and 
odd  subscripts.  The  separation  equations  are  unchanged;  the  characteristic 
equations  of  the  other  three  situations  differ  only  by  the  appearance  of  y 
in  place  of  e  and/or  the  cotangent  function  in  place  of  the  tangent  func- 
tion [Harrington,  1961].  These  geometries  could  be  easily  investigated  with 
a  minimum  of  simple  changes  in  the  program. 

The  cases  involving  losses  both  within  and  without  the  slab  may  be  in- 
vestigated with  the  present  program.  Such  cases  may  well  yield  more  possible 
zeros  to  be  examined,  a  and  6  are  no  longer  orthogonal  in  this  case;  the 
angle  between  them  can  never  exceed  tt/2,  however. 

Magnetic  materials  are  extremely  useful  in  the  microwave  region  and  are 
usually  rather  lossy.  Thus,  magnetic  cases  may  well  be  of  interest.  Unfor- 
tunately, most  such  applications  capitalize  upon  the  anisotropy  of  such  ma- 
terials when  subjected  to  a  dc  bias,  so  that  the  permeability  is  a  tensor. 
A  tensor  permeability  cannot  be  used  with  the  present  program. 

Cylindrical  geometries  may  also  be  of  interest.  It  has  been  proposed 
that  a  lossy  dielectric  coating  be  used  upon  power  lines  to  eliminate  the 
propagation  of  RFI  by  surface  waves.  If  the  results  of  the  present  research 
are  indicative,  such  a  scheme  is  not  practical. 

To  the  best  of  the  author's  knowledge,  there  has  been  no  experimental 
verification  of  either  the  transition  phenomenon  or  of  the  existence  of  the 
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TM,  extraordinary  mode.  The  technology  of  constructing  such  lossy  materials 
is  available  or  soon  will  be;  anechoic  materials  are  commercially  available, 
intended  for  uses  such  as  the  elimination  of  undesired  reflections  in  micro- 
wave measurements.  It  should  thus  be  practical  to  perform  such  experimental 
research. 
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APPENDIX  B 


PROGRAM  LISTING 


The  "$JOB"  card  at  the  beginning  of  the  listing  and  the  statement  numbers  at 
the  extreme  left  of  some  lines  are  included  by  the  Waterloo  FORTRAN  (WATFOR) 
compiler  and  are  not  part  of  the  FORTRAN  card  deck. 
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:JtF  JKc.WlN  =  OtCK  .TIM  =?fPAGHf  =  *lC 

PWOCf".  Af*      PUNI'I    Si.  :      LF    It-  HF   IN AT  It  k.     C  F       It"-      t'wfC«(.AT  ICN     CrfsSIANr*"     i.F 

ILI-FACl       ANC     L ».  4  "  Y      MKC      T-iANSVrPSt      *  *  CNE  T  1  C      'LI'0 
LVlK     A     GHLLMI'-L     ClLL.-CTF.IC     SLAb      INFINITFLY     LCNG, 
INFINITELY     VICL.      ANC      I     C  C  N  T  I  *  b  TC  Fi  S     Tt-ICK,     WH'PF 
f?CTF      TFl      CIM.KHIC     SLAe     ANC     SUF.MCLNCING     M-CIL!' 
.VAY     FAVF     i.lTHt-tv     CLfcCTiSIC     O     VAGNFTIC     LOSSF',     c- 
LI  It-. 
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(.FCICN      I  :      SLF-l-LLNDING     M.CILK.     HAVING     F*  F  L  A  T  1  V  C 
PI'^MITIVIIY       .~1      ANC     RELATIVE     Pfct- 

1'i-AtlLITY     F'L-'  1  . 


X=T 


X^t 


"IICICN     ?:     DIELECTRIC      5L*h     HAVING      SFLATIVi       FFnMTTIV- 
ITY     LH2      AM.     t-UM  IVf.      PL'V^FAPILITY      KU"?. 

////////////////////////////////////////////////////// 
FERFFCI  CLNCLCTCF; 


F  IOU«r   INDICATING   TFF  tFTfiTCY  CF   TFt  RLFFACE  WAVEGUIDE. 

T  I-  F  C-»F.TtSIAN  CCCRCINATI-  A  X H.  S  APE  C  W  I E  N  T  E  C  AS  H  C  L  L  C  W  S  :  THI" 
2-AXIS  IS  CInECTEC  TC  THF   UGH;  TFC  Y-AXIS  IS  DIRECTED  GUT 

TF  THE  PAGE;   THF   X-AXIS   IS  DIFECTCC  I.  P  W  A  F.  C  ,  WITH  X  =  U  AT  THC 
INTF-MFACI-  F'LTViELN  THL  fj  I  E  L^  C  T  R  I  C  SLAb  ANC  FEF-F^CT  CL'i\  CLC  TG>3  . 
TFT  STRUCTLi't   IS   INFINITE   IN  rXT!?NT   IN  FCTF  TF-C  Y  ANC  I 
L  It-FCT  ILNC  ;   TFr   SLAt:   IS  T  CHNTIV..Tf*c   TF1CK. 


fELC*   TI-E   '•LAFj-CCNCLCTtw  INTFHFACt  THF  FIELDS  AhlF  ZFFiC. 

IN  H  e  C I  UN   ?   (   C  <=   X   <=   T   )   THF   F  I  •_  L  C  S  Am; 

H  (  Y  )  =  -Ci»L*CLS(LX)*FXP(-J*FW*/') 

E(X)  =  -Ci  /(  CfECA*t<;  )  *U*KZ«CCS  (  LX  )»FXP  (-J*KZ*Z  ) 

E(2)   -  C2/(J*C^FGA*F;)*G»*2*SIN(UX)*bXP(-J*KZ*?) 

IN     HEGIUN      I      (      T  <=      X      )      THF:     FlcLUS      Ahj 

F  (  Y  >     =  Cl*V*LXP(-V*X)*fXr-(-J*K7*Z) 

E(Xj      =  CI/(L/'t  (..MK|)*V«Ki?*tXP(-V*X)*FXP(-J«K/*Z) 

F.  (  Z  )     =  -Cl/(J*CVLG«*ri)*v«*;*bXF(-V*X)*FXP(-J*Kz*7) 


IN      T»  •=      AHUVt-l     LX*>RFSS  1  t  NS  . 

MY)      IS      THr-     V-CIWECItU     CCPFLNENT     LT      TFT     VTCTGP     H. 

FIX)     ANC     E  <  Z  )      »HH,     NH SPLC T I VFL Y  .      TF'fc     X-     ANC     Z-CIF.GCTCC     CCNTCNfcNTS 
CF     THC     VhCICR     La 
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J      IS      IFF      IPACINAHY     LM  T  :      J     =     CSCR T ( -  I  •  C  )  . 

cmeca    is    r*c    tiiy-s   pi    times   Fhj.iCLCNCY,    rn"    "angular   frequency". 
ti   *md  c?   Afji;    ih-_-   pknmttivitifs  o   thc    t»l   rf.gicns   efecikiuj. 

CI  Af»c  C2  AWF  GCWFLLX  CCl'.FF  IClCNTS  ohlCh  OC  MT  CEP-ENC  UN  X.   V, 
CH  Z,   etl  DL  UFF^NG  CN  TO/Co*.   :  1  .   -2.  LTC. 

U   1=  TFE  X-CCPPCMEIST  CT   THE  FHCPAGATICN  CONSTANT   IN  RFGICN  2. 

J*V   IS   TFT   X-CCCPLMM  LF  IFF   PHCFAOATILN  CONSTANT   IN  RiClCN   I. 

Kl      IS   THF   ^-CC^PC^^^T  CF   TFT  PHCFAGATICN  CCNSTANT   IN  ECTF  SECICNS. 

TFIS  PROGRAM  GCIFRMNES  L,  V.  ANIC  KZ  F  C  H  A  GIVEN  SET  CF  PAhA- 
"FTEF-S.    TFT  VALUES  CF  Z-DIHLCTICN  «  T  I  ^NL  *  T  I  L  N  ,   X-CIRECTICN 
ATTENUATION   (IN  RfGICN  1).  L.  V.   «NC  K  7 .   ALL  EXFRESS.EC  PEw  Eh>EK 
SPACE  HAVELFNGTH,  PLUS  THE  K»TIC  CF  GUICF  WAVELENGTH  TC  FRR;  SPACE 
kAVELENGTh.  ARE  PHINTEC  CUT  ALONG  »ITF  TFF  PARAMETERS  GF  TFI-  CAGF. 

....HMttMII<MI«(MtM.MMt(MltttHt(HHM(ttltl(MH(MttM 

CAfA     WtQUIRFC      (NCTt     THAT     CEFALLT     VALUES     ARE     FHtSENT     FCR     ALL 
FAfcfAyFTF.PS  )  : 

*ITF  TFE  CXCFPTICN  CF   JQATA,  ALL  INPUT  CAflC'  CCNTAIN  A  CHARACTER 
STMING   TITLE   IN  CCLUMNS  2-9,  ANG  FLCATING  FCINT  NUNErRS   IN 
CCLUMNS  20-29,  AC-'*1;,  ANG  60-69',  SCV-  CK  ALL  CF  TFE  NUMFrirBS 
*«Y  HE  ZEHC.   ANG  TITLES  MUST  Eu  LEFT   JUSTIFIEG. 

TFC  FCSSIULE   INPUT  CARCS  AWE : 

CCLUNNS  2-9  2C-i:9  40-A3  60-6<J 

FRGC         (FRGCU^NCY  IN  GF/.) 

SIZE         (THICKNESS   IN  CV.) 

MODCHNGE    (LCHFH  fGCE  NO.)       (UFPER  MCDt  NC  .  )  *SEr  NOTE   I 


ER1P 

PP1PP 
EHiJP 

tm?pp 

KUH  1H 

MUT2P 

VUW2FP 

HAT  IC 

ZST ART 

ID4TA 

DCLTA 

CELUC 

TAriLF 

GLCCTH  IC 

MAGNETIC 

NCELECTH 

NO" AGNF T 

HUN 


(WEAL(FR I ) ) 
(-AIWAGIEH | ) ) 
(HCAL(EH2) ) 

(-A  IfAG(  >-H?)  ) 

(  REAL  (MUH  1  )  ) 

(-AIMAGIMURI ) ) 

(HEAL (NCR  I  )  ) 

(-A  I*A-,(*LH?)  ) 

(  T/LAMiDA_C  ) 

(RCAL(/))  (MKAC(Z))     *SEE     NOTE     ? 

(NL.     GF     VALUES)      *S     G     NGTF     3 

(GCLTA) 
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STUP 

NOTE   i:   Tllf~SF   VALLlS  Ak|_  CONVERTED   T  (.   IMFll  PS--*  FLCAT1NG  PlIlNT 
NUNUER  SL1CHTLY  LARGER  THAN  THE   INTCGCR  SHCULC  Hk     O'iFD  TU 
KNSUHf   THAT   THE  PRCPFR   IMKiH  KILL   bF  UHTAINLD.    THL 
INTFCERS  MUST  IE  EVEN. 

NOTE  2t  Z  HIRE   IS   Tt- ':  FWCGRAW  VARIAEL-  U*T.  NCT  THE  CCCRCINATf- 

VALUE.    RLAL(ZSTARI)   SHULLC  !  ".   >=   C   TC  AVCID   AN  UNLNMNG 
LCCF  ft  I  TH  IN  HAN. 

NCTF  2:      THIS  VALLf   IS  ALSf.  CfNVLt-lci;  TC  AN   INTEGER.  SC  THE  PRE- 
CAUTION Lf  NCTE   1   AP,~LItc.    T  H  F   INTEGER  NAY  NCT  fair  LA^r.ru 
THAN  7S  UNLESS  iDATA  IS  DIKtNSICMrC  ACCORDINGLY.    IMWFCIA- 
TF.LY  FCLLLIAlNC  THIS  C*n;:  wucT  Hit  A  CARD  CR  CAMUS  CCNTAIN- 
ING   THC  VALLCS  LF   ICATA   IN   10Fr:.0  FCRKAT. 

EITHFR  FREC  AND  SIZE  LR  RATIC  WAY  cF  USED.    IF  FREC  AND  S  I  Z  fT  ARC 
PROVIDED.  HAT   IS  CALCULATLC  HRC*  Tt-rf. 

EW1F  1H'<L  NCRi-PP  SET   THt   VALLF  CF   THAT  RARAKFTtfR. 

DELTA  SETS  THt   INITIAL   INCREMENT  VALUE  FLR  STEEPEST  CESCENT. 

*CATA  DETERMINES  THE   VALLtS  CF  L\fl2FF  (R  NLfl^PF  DURING  PASAyfllE" 
SteFTPS. 

ZSTART  PROVIDES  A  STARTING  VALUE  CF  THE   VARIAELE   ITERATED  UPON 

FUR  THT  FIRST  PASS  *H  I  CH  K«Y  be  USED  TC  CVERRIUE   THt  VALUE 
CmLCLLATFO  HITHIN  THF  PRCCWAW   IH  NhCSSSAPY. 

i  lectr  ic/ncllfctr ic  ccntrcls  hhi-then  lr  nlt  elfctric  lcsses  are 
Swept  . 

magnf t  ic/numagnet ic  ccntrcls  kmtthlfi  cp  nct  magnetic  lcssfs  arc 

swfpt.   if  hlth  typts  of  s  *  t  >r  h  are  specifiec  the  flcctpic 
s  m  *=  e  f  is  performed  first. 

talile  prcdlces  tahular  cuthlt;  dellg  produces  an  cutput  *itf 

thc  derivatives.  efihcr  values,  anc  z  values  tc  assist  in 
debugging. 

upon  cncountfminc  the  run  card  tht  p«cg«ah  is  executed  using  thf 
oaramet'rs  in  effect  at  that  tine. 


UPON  ENCOUNTERING  THE  STCP  CARD  SLHRLUTINE  EXIT  IS  EXECUTED. 

UPON  cNCOUN Tl  RING  ANY  TITLE  NCT  L  I  S T L L   IN  THF  ARRAY  LALEL.   THE 
LNOtCUOEAeLF   TITLE   IS  FwINTFC  CUT  ANC  EXECUTICN   IC 
AHANDLNED . 

MMIM>MMIMlM<MHIMI.<<M(MHtl<<M>t.tt (...Ml. 


CCMFLCX  Z.  ER.  K1SC.  K2SC,   J,  L,   V,  C  S  C  -i  T  ,  KZ,  ZTAN 

COMFLEX  DFLF,  ZC 

COMPLEX  ZCK,  STCRE 

COMFLFX  ZSTAHT 

COMFLFX  OZ(?A) 

REAL     MUR1P,     MUR1PP.     MLR<>P.     MLR2PF 
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CblALLT     PAwAM.  TFPS 

CATA     FCH7/3.C/.      T/l.C/t     NL/l/t     NF/1/.     F  W  I  F  /  1  .  C  /  .     LH1PP/O.C/, 
&  FF2P/2.C/,    '.3U2FF/C.C/.     MLR1P/1.C/,     NLi".  1  FF/C  .0/ ,     PlRcP/l.C/« 

&  M  i-U'PP/t  .  r  / 

HATA      tCATA     /f.C.      C.«r.     C.3.     0.'.,      C  .  ..  ,     C.6.     (..7,     0  .  P  .     0.4,      I.e.      I  .  «•' • 
C         1.4.      I.e.      l.P,     2.C.     ?.j.     3.C.      3.ti»     4.C.     'i.l:,     5.0«     5.5,     £.0, 
C        bi  *C.C  / 

C  EF  T  N  t     FLNCTlt.N     IT '1  HAT  Ft;     I.  P  C  N 

F  </  )     =     CAbC(Z-f3H*RAT*TfcCFI*CrCHT(  (Kir.C-Kl?C)/(C«**2+ZTAN(/)**r  )  )  ) 

IMTIALI/F     CCNSTANTS     AMJ     SfclTOtrS 

I  F  E  *>     =     0 

i  f  r  *    -   o 
I  7  S  »    =    f ; 

NPASS     =     0 

IFCR*     =      1 

NCT     =     23 

J=( C  .  .      1.  ) 

TViuF  I  -<  .  *3  .  14  lt,<527 

PI  -J .  14  1  J9^7 

IrLl'C     -     CCS*PI 

«AT     =     ICHZ*T/^f  .C7<ij 

»n  I  T  F      (3,     ?4<5C  ) 

teHITF      (3,     <J4<;  1  ) 

* ill  T  F      (  3  .     9  4  5  2) 

*PlTt      (3.     <;4<iJ) 

»HITF      (3,     FCC2) 

IF     TFif     PWCGfAV     IE     LS^r.      IN     FCSlniAN      IV,      n,     CC     C^C     MJST     bfc     SU'PLI'iC 
€  FC!!     GMT      'I 

FCC     HEAT.      (1,     PCCC.     l;NC-(U)      (CHAH(I),      I      =      1,     2  0) 

*«1TF.      (4,      8CCC)      (OAR(I),      I      r      l,     2  c  ) 

*MTF      (3,     L'CCl)      (CMR(l)i      I     =     1.     2C) 

CC     TC     HOC 
6C 1     WE* INC     4 
CCC     WEAL      (4,     <3C<5<:)      [NVAWI(l),      |NVAKI(2),     VALUTA,     VALUtl.     VALGFC 

KJs'F      Tf  5TS     MAY     dE     ACLEC      IF      TFF     CATA     STATizCn'NT     LIST      IS     ADCtC     TG 
t  AND     TFE     CIPLNSICN     ARGLCcNT     CF     LAflCL      IS      INCREASEC     ACCCRCINGLV 

IF      (       INVArtl(I)      >LC.     LAMcL(Cl)      .ANC.      I N  V  A  K,  I  (  2  )      .EC.     LAefcL(C2)      ) 
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4C 

4  1 

42 

4j 

44 
4E 
4t 
47 
4f 
4? 
EC 
E  I 


£4 
c  c 

Et 
£7 
E£ 
=  c 

er 

6  1 


I  F   (   IKV*^  I 


I F   (   INVAR  I 


IF   (   INVAH [ 


IF   (   INVARI 


IF   (   INVAR  I 


IF   (   INVAR  I 


IF   (   INVA<I 


IF   (   INVAI'l 


IF   (   INVARI 


IF   (   I  N  V  A  M  I 


IF   (   INVARI 


IF   (   INVARI 


IF   (   INVAW I 


IF   (   INVAR  I 


IF   (   INVARI 


IF   (   I  N  V  A  W 


IF   (   I  N  V  A  .-<  I 


IF   (   INVARI 


IF   (   INVAR  I 


IF   <   INVARI 


IF   (   INVARI 


IF   (   INVARI 


IF  LAP^L   IS  NCT 
HfilTF  (3,  ooof  ) 

»h  i  te    (j.    <;cci) 

*W  I  IE  (3,  <5Cf  2  ) 

mnr  (3,  octi) 

*RI  TF   (3.  0CC4  ) 
STOP 
0  C  1  FCh^  =  V  A  L  C  E  A 

RAT  =  FC,t-/*T/20 

cc  ir  occ 

0C2  T  =  VALUEA 

GAT  =  FGI-^*T/2  0 
GO  TC  OCC 


CC  TG  '.CI 

(  1  )   .  1:  C  . 

L»RHL( C  3 ) 

.AM. 

I N V  Ak  I 

(2  ) 

.EC. 

LAPEL <C4J   ) 

GC  IC  -»G? 

(  1  >   >f-Gi 

LA(JtL(C5) 

.ANT.. 

I  N  V  A-.  I 

(2) 

.EC  . 

LALtL(Ct)   ) 

GC  TG  «i 

(  1  )  .EG. 

LMlL  (C7  ) 

.  ANC. 

INVAh I 

(2) 

•  EC. 

LAFdL(Cr)  ) 
G  C  T  C  «i  C  4 

(  1  )  .FC. 

LAILL(C) 

.  ANC  . 

INVARI 

(2  ) 

.EC. 

LAFkLdO)  ) 
GG  TC  '.C  • 

(  I  )  .KG. 

L  At'FL  (11) 

.  ANC  . 

I NV AG  I 

(2  ) 

.EC  . 

LAFcL (12)  ) 
G  C  T  C  H  C  t, 

(  1  )  .KG. 

L AP^L (13) 

•  ANC. 

IKVAmI 

(2  ) 

.tC. 

LAPJL (  14  )   ) 

GC  TC  007 

( 1 )   .  '  C. 

LAI-LLd'i) 

.  ANC  . 

INVAh  I 

(2  ) 

■  EC. 

L  Ah>L  (U)  ) 
GC  TO  SCC 

( 1 )  .EC. 

LAHH.I  17  ) 

.  ANC  . 

INVARI 

(?) 

.EG. 

LAe;L( lb)  ) 
GC  TC  «'< 

(  1  )  .  L(.  . 

LALtLI  1  0  ) 

•  ANC. 

INVARI 

(  2  ) 

.EC. 

LAPEL (2C)  ) 
G  C  TO  0  10 

(  1  )  .EC. 

L  Af  tL  (  2  1  ) 

•  ANC. 

I  NVAl.  I 

(2  ) 

•  EC. 

LA EEL (22)  ) 
CC  TC  Oil 

(  1  )  .  K  C  • 

L  AhtL (23) 

.ANC. 

1  NV  A  ■«  I 

(2  ) 

.EC. 

LARr.L(24)  ) 
GC  TU  012 

(  1  )  .EC. 

LAPEL  (25  ) 

.  ANC  . 

INVAh  I 

(2) 

•  EC. 

LAEtL  <2t.)  ) 
GC  TG  il.i 

(  i  )  .  r  c . 

LAfc'GL  (  27  ) 

.  ANC  . 

INVARI 

(2  ) 

.EC.  . 

LAEEL (26)   ) 

CC   TO   0*4 

(  1  )   .LC. 

LAf'EL  (20  ) 

.ANC  . 

INV  Ao  I 

(2) 

.EC  . 

LAEEL  (JO   ) 

GC  TG  014 

(  1  )  .EC. 

LAEEL (31 ) 

.ANC. 

I NV AK  1 

<?  ) 

•  EG. 

LAFEL(32)  ) 
G  C  T  C  OIL> 

(  1  )  .EC. 

LAEtL (33 ) 

.  ANC  . 

INVAR  I 

(  2  ) 

.FC  . 

LACEL(34)  ) 
G  C  T  0  0  I  b 

( 1 )  .tc. 

LAfiFL  (3=  ) 

.ANC. 

IN  V  Ah  I 

(2  ) 

.EC. 

LAfccL(Jh)  ) 
CC  TC  R17 

(  I  )  .CC. 

LAEtL (37 ) 

.  ANC  . 

INVAR  I 

(2  ) 

•  EG. 

LAPEL(JP)   ) 

G  C  T  L  411 

(  1  )   .EC. 

LAC'EL  (  30  ) 

.  ANC  . 

INVAR  1 

(2) 

.EC. 

LAEtL  (HI  ) 
GC  TG  or; 

( 1 )  .EC. 

LAPLL (4  1) 

.ANC. 

1 N V  Ar  I 

(2) 

.EC. 

LAEFL('i2>  ) 
GC  TO  OiO 

(  1  )  .EC. 

LAEEL  (43  ) 

.ANC  . 

INVAR  I 

(2  ) 

•  EC. 

LAECL (44  )  ) 
GC  TG  0^1 

(  1  )  .EC. 

LAPEL  (4'o  ) 

.ANC. 

INVARI 

(2  ) 

.EC  . 

LAEEL  ( 46  )   ) 

GG  TG  0,  > 

IN  LIBRARY,   V.RITE   A  VtSfAG'. 
INVARI(l).   1NVA«I(2) 


ANC  STGF 


0  7  0  3 


0  70. 
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74 
75 
7t 
1  7 
76 
7  5 

ec 

6  1 
f  2 
C2 
£4 
tl 

e< 

£7 
€£ 
£5 
9C 
5  1 
52 

52 
54 
5e 
5t 
57 

5£ 
55 
ICC 
1C  I 
1CS 
1C2 
1C4 

ic; 
ice 

IC7 
ICC 
ICC 
1  IC 
1  1  I 
112 

1 12 

1  M 

1  is 

1  It 

I  17 
1  16 
1  15 
ISC 


12  I 
122 

122 
124 


IS! 
ISC 

127 


5C2 

SC4 
5C5 
5Ct 
5C7 
9Ct 
5C5 
910 
511 
9  12 
512 
914 

51  = 

5  16 
917 
5ie 
515 

52C 

52  1 

522 

524 


NL  = 

NH  = 
CL  T 
It  S* 
CO  T 
ICS* 
GC  1 
FM  IP 
GC  I 
FB  ll- 
CC  T 
F«?P 
GU  1 
fcH?P 
CL  T 
VUW1 
GO  T 

cum 

GC  T 
NCK2 
CO  T 
MJW2 
GO  1 
RAT 
CO  T 
t?  14 
IZSin 
GO  T 
IF  ON 
GO  1 
IFDH 
GO  T 
ItS* 
GO  T 
I  VS» 
GG  T 
CALL 
NCI 
HEAC 
GC  T 
CELT 
GO  T 
trt=( 
K  ISC 
K2SC 
CD  2 
I»  = 


VALLFA     4      l.COl 

VALUTI       ♦      l.CCl 
C     5CC 

=      1 
C     5C0 

=     I 
C     5CC 

=     VALLt"» 
C     <5CC 

F     =     VALLFA 
C     5CC 

=     VALLTA 
C      5  CC 
F     =     VALLrA 

t    <icr 

f  =     VALLEA 

L  500 

FP     =     VALLEA 

C  5  CC 

F  -     VALLFA 

C  r,CC 

pp     =     VALLEA 

C    ^cr 

-     VALLFA 

r    sec 

R!     =     VALLI.A     ♦     J^VALLO 


c  5 

fc  - 

c  5 

*  = 

C  f; 

C  <> 


Ur.A 

9C57)      (tCATA(I),      I      -      1.     KCT) 


C      <i 

EX 
-     V 

(4 
C     5 
AC     =     VALLFA 
C     9 
FHi 


ICC 


ICO 

>p- 

(  N 


j*riJt)pp)/(i.HiF-j*rwipp) 

UHlP-J*fLH IPP)*(£H1P- J*e« IFF ) 

UWPP-J*WL^2PP  )*(EH2  P-J*Fhi?FP) 

M.     INh  •     2 

I 


5£C 


CHCCCE      STAHTINC     VALLT      FCH     / 

IF      (      US*     .FC.     C      )         CC     TC     SSC 
I     =     i>STAHT 
IZ?.*     =     C 
GO     TC     92c. 

TEST     Ft.l?     CLTCFF 

HACILS     =      f  V»CF  I  *HA  T+SCPT  {  KLR2P  »F  W2  P- f*L^  1  F  ♦  E  P  IF  ) 
IF      (     RADIUS     .GE.      IHFI/2.      )         GC     TC     55b 
/.     -     I  IH     -      1     +     C  .  1  )  +  P  I  /  2  .      ♦     J  ♦  0  .  j 
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ISC 

12S 
12C 
121 
122 
122 
12« 
128 
12< 


127 
12£ 
1  2< 
I4C 
14  I 
142 
142 


1  44 
14S 
14£ 
147 
14£ 


14C 
1£C 


1£  I 
1S2 
1£2 
If  4 
Iff 

l£t 

if? 
i?t 
if? 

itc 

16  1 
If  2 

it:- 


1  t  4 

let 

1C7 


ICt 


GO     TC     S?£ 
Sff      IF      (      t)ACILc  .Gl.  .      (  I*+  1)*PI/P*       )  jt       T<       9t>C 

Z     =     C«*f4>MAnit£     ♦     C  .  1*  Ix*Pl /?. 

CO    IT     t2') 
SCC    Z     =     C  I»*,+  0«9)*HI/2!« 
S2=    KPA55     =     1 

uxrc    =   c 

1  mac    =    c 

IF     (      li  r\\     .  i  r.  .     0     )        GC     IC     9  3  C 

?»ITF     I  LECTH  IC     LCS  £Lf 

Sit     lfrLCC     =      IFLt-C     ♦      1 

FHPFF     =      TDAT A<  If  L  fc  C  ) 

KicC      -      (MjH,F-J«fth<:FP)Md^P-J*tf.;,FiM 
fcW=(FRPF>-J4FF<?PF)/(trniF-J«EF;iFP> 
GO     Tf.     <>3=i 
S2C     IF     (      If<S*     .Kt.     C      )        CC     TC     SJb 

IF      (      Nt'A";     .Cf.      1      )         2     =     7JIAWT 

r  Vk  i ;  hF     f>Ar,NFTIC     LCSif.  2 

S31      IK  AC-     =      If*  AG     ♦      1 

PU'JSFF     =     tHATA(lfvAC) 

K23C     =     (KLH2F—  J*H»LB2PP)*(t:K2:P—  J*PK2FP) 

S2£    COMING- 
IT     (      !HJ>»f*     .cU     C      )        URITE     (3,     SSfcl) 

ITEPMION     FFCINS     f-Ff<t 

DC      I     K     =     1  ,      ICt 
FLAG     =     C 

T«Y     LINf  »^      ITLkATILN     AKU     CHECK     fcFCTFFP     C  K     NLT      IT      IS     USEFUL 

ZCK     =     EP*RAT*TV»(JFI*CSCHT<  (K2S0-KISC  )/<Fh!**i.*ZTAMZ)**2)  ) 
IF      (       IERP      .FG.      1      .AND.      IFCWN     .  F  G  .     C      )         **ITF      (3,     9401) 
FtMf.     =     K  <  /  ) 
FCKFKFS     =     f(/CK) 

IF      (     rCKFWtv     .f,T.     Ft-Ffc         .CW.         CAeS(Z-ZCK)      .GT.     O.C3*CAES(Z)      ) 
G  GC      IC     7 

STO«F      =     ZCK 
ZCK     =     / 

2  -     STCWF 
FFHH     =     F  (  Z  ) 
FCKEFB    =    F(7CK> 

rLAc    =    i 

7  IF   (  FFPW   .LT.   l.t-C7*BAT   .AMI.   [FLFf   .  F  G  .  C   )    *RITE   (3.  ',<»  -C  ) 
IF  (  FfcHG  .LT.  l.E-C7*RAl   )    GT  TC  4 

GET  CewlVATIVtC  FCf*  £TEtf-E'ET  DESCENT 

FPX  =  (F(,-  +  C.CCC  1  )-F  (  Z-U.CCO  1  )  ) /C.tCC? 
FPY  =  <F<Z+J*C.CCC  1  )-F < Z-J*C .CCC 1 ) )/G.CGC2 
CELF  =  FPX  +  J4FFV 
t.HLffC     =  CAHS(CELF) 

IF  CFLMAG   IS  ZFPL.  PGSt*T   TC  THE  CAPtlACi  CAN  LID  AFFBCACh 

IF  (  OFLKAG   .EC.  C.   )    GC  TC  b 
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1  7C 

n  i 

17J 
173 
1  74 
17? 


17< 
177 
176 
1  ~i< 
ICC 

1(1 

IE2 


162 
164 
16* 

let 

167 

ice 
it? 

i<;c 
l?  i 

ISi 

is: 
i«« 

1SE 

is< 

1<!7 

i«e 

l?c. 

2CC 
2C  I 
2C2 

2c: 

2C4 

2ce 
2ce 

2C7 

2ce 

2CS 
2  K 
2  1  I 
2  l< 
213 
214 
21S 
?  \t 
2  17 


2  It 


RZ     =     Rl  AL  (  Z  ) 

A  IZ     =     AINfiCt/) 

IF      (      lfV<t       .t'C.     C      )         »WMt      (3.     S<.t^)     FFHh,     FLAG.     FCKEfH.     rt2«      A[/, 
t  F  f-  /  .     F  r  Y 

IF     (     C  F  L  V  A  G     .<jT.      l.C      )        CtLF     a     CFLF/ CELIAC 
OrLTA     =     DEL T AC 
KC     =     C 
20     =     Z 

THY     5THPFH      LF'JCt-NT 

t     Z     -     ZC     -     OCLTAtDFLF 

IF      (     F<7)      .LI.     TFI.I)      )         oC      TC      1 

CELT*     -     Ot:LTA/i. 

KC     =     KC     ♦      1 

IF      (     K  C      .LI.     2  S      )         G II      T  C     c 

7     =     70 

FLAG     =     FLAG      ♦     2 

FCGIK     GARhAGF     CAN     LIC     ATTFKPT 

6     FLAC     =     FLAG      ♦     d 

KC2     =      C 

HFLZ     =     C.CCO I 

IF      (FFHH      .LT.     IrLZ      )         CtLZ      =     FrJfi 
2  1  C     C  C     2  (  C      1 1 1 Z     =      I  .     ?  4 

TFt-IA     =     T*CP  l*(  ICZ-  1  )  /?4 

CZ(irZ)     =     CFLZ*  (CCS  (  TFFTA  )♦  J*S  IN(  THE  TA  )  ) 

£HR(1U/)     =     F(/+DZ(  ICZ)  ) 
3CC     CUNT INLh 

10/     =     1 

IMIN      =      I 

F)nii  in    =    ru -I  (  I  ) 
2C  I      IDZ     -      ID/     ♦      I 

IF      (      lrR«(ICZJ      .LT.      FHHMK      )  IWIN     =      IC7 

IF      (      ^RH(IOZ)      .LT.     FNRMN      )         ERR*  I  It     =     »'"^I  IMM 

IF     (      IUZ     .LT .     24      )         GC     IC     3CI 

F  IMIN     =      IM  IN 

IF      (      HRRM I N     .LT.     fr.Kti      >         GC     TC      302 

FJFfLZ     =     DFLZ/2 

KC2  =  KC?  ♦  I 

FKC?  =  KC2 

IF  (  KC2  .LF .   1 7      )         GC  TC  310 

dZ  =  RFAL(Z) 

A  1 Z  =  H  I W  A  G  (  Z  ) 

IF  {   IFU'-M  .tG.   1   )    GC   TC  4 

VI9ITF      (3.      59>:2)     Ftl>ii.     FLAG,     FCK.rKR,     WZ.     AlZ,     FIMN,     FKC2 

WR  I  TC      (3.     99^1) 

GG     TC     4 
2C2     I     -     2     ♦     DZIININ) 

IF      (      [  fc  RR      ,r-C.      1      .  A  N  0  .      IFCRS      .l.C.      C      )         WhlTC      (2.     S  4  C  1  ) 

RZ     =     WFAL(Z) 

A  1/     =     A  IN»AG  (  Z  ) 

IF      (      IFGHf      .FO.      1      )         GC     TC      I 

\»«ITfc      (3.     <;<3£2>     Ft>S.     FLAG.     FCK~fc^.     HZ.     AI7.     FIWIN.      FKC2 
1     CCMIMf 

END     CF      ITcHATICN     LCCP 

»H  I  TF     (3.     <;<5eo) 
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2  IS 
22C 
21  I 
2ii 
222 
224 
225 

sit 

221 
22t 

2?  = 
2  2C 

21  I 


222 

222 

?:-4 

225 
22f 

2  2  7 
22t 
22? 
24C 

24  I 
245 
2  4  2 
2*4 
245 
24t 
24? 
24F 
24C 

25C 

25  1 
25J 

252 
25  4 

25t 
254 
25  7 
25fc 
25  c 
2CC 
2t  I 
2ti 
?c2 
2'.  4 
2f  5 
2*  t 
2tl 
2tt 
2€« 
2  7C 
2  7  I 
27! 
272 


CUM  I 
L  -  I 
M  -  i 
K  J  i 
LJ  AT 
ATTX  = 
A  I  I  l- 
(j^  =  H»: 

te^A  i 

L  LJ  =  A  | 
L9  =  Gt 
C  IC  =  A 


KLf 
/.<AT 
♦/ T AN ( 
CC.G«T  ( 
=     Wt  AL 

-fi.  t  (    5 

AL  (C  ) 

NAO(l>) 
AL  (  V  ) 
MC,(V  ) 
AL(K/) 

I  ,v  a  r  (  x 


l  )  /  ( ■_  r-  ♦  p  a  T  ) 

A.O*P[«PIM»l.MIF-JtMLK  IFF)*(FKilF-J»EHlFP) 

(k/  >/r»cf i 

F-  0  c  3  ' '  ♦  F.  -"  A  L  (  V  ) 
tfr<je3c*Al*AG(KZ) 


V»*?  ) 


/ ) 


cnnr?'"    cltflt    rcKNAT    HfKt; 


IF      ( 

X  =  <1  (:  A 
Y  =  A  I* 
IN  HI  T  r 
»,!  I  If 
XCK      - 
VCK     = 
»RITF 
»>-'  I  Th 
IF      ( 
h  A  =  t*F 
F  F>  =  A  I 
FC      = 
FC      = 
**  I  TF 

»uiir 

►  R  I  T! 

fcRITE 

*i<  I  TF 
CL      Tf 
4CC      IF      ( 
»R  I  If 


IF     ( 
Vt*  I  Tt 
CL      TC 
IF      I 
IF      ( 
teW  IT  t 
kRI  Tfc 
CC      If 
IF      ( 

•  •Mir 

IF      ( 

NP«?5 
IF      ( 


IFC  ■> 
L(  /  ) 

AC.  (  I 
(  3. 
(  3, 
RfcA 
A  I  V 
(  3. 
(  3, 

ic^h 

AL  (  i 

yA(i( 

Rt  AL 

A  I  V  A 

(    3  . 

(   J  . 

(    3, 

(  3. 
P.      R 

(  3. 

4  5C 
M'/\c 

(  3. 

(  3. 


) 
<5  <5  «i <i  ) 

L ( /CK ) 

AC  (  /.C« 

<:.  c  F  4  ) 

.ir,, 

*/TAM( 
/ */T AN 

( •-:  r  *  c  s 

b( p^ *c 

<;<;<;<•') 
i-v;(: ) 
<;  <;  '■■  b  ) 

AT,      LI! 

<=<;<;  4  ) 
•     .c  i . 

94Cq  ) 


1      )         GC      TC     4CU 


K  »     KC  .     OtLT A 


ir-;,     [;■   LI'AC 


XCK  ,      YCK 
I       .AM..      [FLh'V      .  12  C  >      C      )  *MTr       (3,      •'3401) 

/  )       ) 
(  /  )  ) 

CRT(  (T*CPI*RAT)**2*IKiJEG-KlSC>-Z**2)  ) 
<"C^T((Tv«CPl*r<AT)**c'*(i;'"<.-Ki'G)-^**i.')) 

FA,  r t .  r c ,  rn 

Fb'KK  ,  FCKFC^ 
[  F» 

l:R1F.  cklFP,  PIR1F,  M»IFP,  FR2P.  IrHPFP,  NC 
A  I 
C  *_  .  Cfc,  C7,  CP.  C  «i  i  C I  C  ,   ATTX,   ATT/ 


:m2F,  t  w  i.'  F  P  ,   Mj. 


40 
AC' 

4i( 


IF  ( 
IF  ( 
IF  ( 
IF  ( 
CCN1  I 
CO  TC 


II  LLC  .GT. 
(3,  S4r  7 ) 
ASC 

ii  Ltc  .re. 

I^L^C  .GT. 

(i,  «;4<;(. ) 

(  3,   .4'=,  =  ) 
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ABSTRACT 

The  purpose  of  this  research  is  to  examine  the  effect  of  losses  upon 
the  propagation  of  surface  waves  and  leaky  waves  over  an  infinitely  wide, 
planar,  electrically  thick  dielectric  slab  coated  onto  the  surface  of  a 
perfect  conductor. 

A  transcendental  equation  having  complex  coefficients  in  the  general 
case  is  solved  in  order  to  obtain  the  propagation  constants.  A  computer 
program  has  been  written  to  perform  this  task.  This  program  relies  upon 
the  numerical  techniques  of  linear  iteration  and  steepest  descent. 

The  quantities  of  interest  are  the  attenuation  in  directions  parallel 
and  normal  to  the  surface  of  the  slab,  the  guide  wavelength,  and  the  angle 
at  which  the  wave  enters  or  leaves  the  slab.  The  results  of  the  computer 
analysis  indicate  that  for  a  surface  wave  mode  there  is  an  abrupt  transi- 
tion in  the  behavior  of  each  of  these  quantities  occurring  at  a  particular 
slab  thickness.  In  the  attenuation  parallel  to  the  slab,  this  transition 
is  from  decreasing  attenuation  at  high  losses  for  a  thin  slab  to  increasing 
attenuation  at  high  losses  for  a  thick  slab.  Furthermore,  it  is  found  that, 
for  sufficiently  high  losses,  odd-numbered  TM  modes  exhibit  the  behavior  of 
a  surface  wave  though  such  modes  describe  physically  unreal  waves  in  the 
lossless  case. 


